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SECTION  ONE 


INTRODUCTION 


As  a  signal  is  transmitted  from  one  point  to  another  above  a  surface 
with  electrical  properties  differing  from  that  of  the  propagation  medium, 
reflected  signals  may  be  generated.  These  signals  then  appear,  along  with 
the  signal  travelling  over  the  direct  path,  at  the  receiver.  These  additional 
signals  directly  affect  the  communication  and  direction  finding  capabilities  of 
a  surveillance  link  operating  in  this  environment.  These  signals  are  com¬ 
monly  called  multipath  and  it  is  the  purpose  of  this  report  to  describe  models 
of  such  phenomena  that  are  useful  for  ascertaining  the  efficacy  of  such  a  link 
for  surveillance. 

Specifically,  the  models  developed  are  to  be  used  to  evaluate  multi¬ 
path  effects  on  the  DABS  (Discrete  Address  Beacon  System)  performance. 

For  this  system,  a  downlink  receiver  (sensor)  is  cn  the  surface  and  the  trans¬ 
mitter  is  located  in  an  aircraft.  Clearly,  due  to  reciprocity,  the  role  of  the 
transmitter  and  receiver  can  be  interchanged  but,  in  this  report,  we  shall  use 
the  stated  nomenclature.  The  link  being  discussed  is  one  where  the  receiver 
is  close  to  the  ground  (5-50  m  high)  while  the  transmitter  is  at  a  highc-  alti¬ 
tudes  (500-15,  000  m)  and  the  range  is  large  (10- 100  km).  This  is  in  contrast 
to  ground-to-ground  or  air-to-air  links  as  found  in  mobile  radio  or  satellite 
communications  respectively. 

What  we  shall  do  in  this  report  is  consider  those  effects  which  are 
relevant  to  the  ground-to-air  environment  and  develop  models  which  are 
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appropriate.  Our  interest  will  not  be  in  the  microstructure  of  diffraction 
effects  due  to  specific  structures  but  will  be  towards  developing  models  that 
can  be  easily  used  to  evaluate  the  multipath  effects  on  system  performance. 
Namely,  we  shall  be  interested  in  the  general  macrostructure  of  the  multipath 
environment. 

The  two  main  functions  of  this  report  are: 

1.  Model  Description:  In  this  area,  we  evaluate  the  topo¬ 
graphical,  geometrical,  electrical,  and  system  para¬ 
meters  and  how  they  affect  the  RF  signal  at  the  receiv¬ 
ing  antenna.  Specifically,  specular  and  diffuse  multi- 
path  are  discussed  in  detail  in  terms  of  the  above 
parameters. 

2.  System  Performance:  Once  the  multipath  model  has 
been  developed,  the  effects  of  the  multipath  signal  on 
the  SNR  (signal -to-noise  ratio)  and  direction  finding 
capabilities  of  a  monopulse  antenna  with  a  2°  beam- 
width  and  vertical  sector  beam  are  ascertained. 

1.1  REPORT  OUTLINE 

In  Section  2,  we  discuss  first  the  geometry  of  the  aircraft  and  the 
sensor  by  introducing  a  canonical  coordinate  system.  The  large  scale  topo¬ 
graphic  features  such  as  hills,  buildings,  etc.  are  then  defined  in  terms  of 
this  coordinate  system.  We  then  proceed  to  discuss  the  difference  between 
specular  and  diffuse  multipath  and  the  surface  conditions  that  give  rise  to 
each  of  them.  A  model  of  the  RF  envelope  of  the  receiving  antenna  output  is 
presented.  This  equation  is  the  basis  for  the  multipath  model  development. 


Section  3  discusses  the  nature  and  effects  of  specular  multipath  on 
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the  received  signal.  Using  the  geometrical  models  developed  in  Section  2, 
a  geometrical  optics  analysis  of  specular  reflections  is  developed.  The  path 
length,  angle  of  incidence,  azimuth,  and  elevation  of  the  received  specular 
multipath  signal  are  determined.  Since  specular  multipath  arises  from  large 
scale  surface  areas,  the  reflection  coefficient  is  studied  in  detail.  The  effects 
of  electrical  properties  at  the  reflecting  surface,  finite  aperture  diffraction  and 
losses  due  to  diffuse  multipath  are  considered.  The  section  concludes  with  a 
discussion  of  the  use  of  the  models  in  the  analysis  of  signal  fading  and  monopulse 
azimuth  estimation.  These  issues  are  related  to  the  performance  of  DABS. 

The  model  for  diffuse  multipath  is  developed  in  Section  4.  In  contrast 
to  specular  multipath,  which  is  a  highly  deterministic  phenomenon,  diffuse 
multipath  is  random  in  nature.  To  account  for  this,  a  model  is  developed 
for  the  second  order  statistical  properties  of  the  multipath  field.  A  function 
called  the  channel  spread  function  is  introduced.  The  spread  function  is  then 
evaluated  in  terms  of  the  surface  properties  using  the  Kodis-Barrick  scatter¬ 
ing  cross-section  model.  Section  4  concludes  with  a  discussion  and  evaluation 
of  the  effects  of  diffuse  multipath  on  azimuth  estimation  performance.  It  is 
shown  that  there  is  a  negligible  effect  in  the  DABS  case. 

Section  5  discusses  other  effects  that,  although  not  included  in  the 
models,  may  in  time  be  important.  These  effects  include  refraction,  line- 
of-sight  limitations,  shadowing,  sea-surface  scattering,  and  composite 
scattering. 

The  last  major  technical  section  is  Section  6.  Here  we  discuss  the 
various  applications  to  which  the  models  may  be  put  to  use.  Specifically, 
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we  mention  link  and  azimuth  estimation  performance  and  localization  techni¬ 
ques.  Also  discussed  are  antenna  design,  signal  processing,  and  siting. 

The  conclusions  of  the  report  are  detailed  in  Section  7.  Basically, 
we  show  that  DABS  should  function  acceptably  based  on  the  multipath  model 
developed  in  this  report.  We  also  comment  on  the  other  uses  of  the  model 
and  briefly  on  its  shortcomings. 


SEGTION  TWO 


FUNDAMENTAL  PHENOMENOLOGY 


In  order  to  fully  understand  the  effect  of  multipath,  on  the  energy 
received  at  an  antenna  aperture,  it  is  first  necessary  to  describe  the  total 
environment  in  which  the  signal  transmission  is  taking  place.  This  environment 
includes  the  reflection  geometries  involved  between  transmitter  and  re 
the  nature  of  the  waveform  transmitted  and  the  effects  of  the  receive*-  antenna 
on  the  received  signal  energy.  In  this  section  we  shall  treat  each  oi  these 
issues  separately  and  on  a  level  which  will  provide  the  reader  with  a  qualitative 
as  well  as  quantitative  understanding  of  the  problems  involved. 

2. 1  SENSOR  GEOMETRY 

The  basic  geometry  of  the  aircraft  and  DABS  ground  receiver  is 
quite  simple.  In  order  to  avoid  certain  analytical  problems  we  shall  first 
assume  that  we  have  a  flat  earth.  Let  us  then  choose  a  point  O  on  the  earth 
plane  and  at  this  point  construct  a  set  of  three  orthogonal  axes,  one  being  nor¬ 
mal  to  the  plane.  This  will  be  defined  as  the  canonical  coordinate  system  and 
will  be  used  throughout  the  report.  The  interrogator  antenna  DABS  sensor  is 
located  at  an  altitude  z^  above  the  origin.  This  point  represents  the  phase 
or  geometrical  center  of  the  antenna.  It  is  that  point  to  which  all  phase  delays 
are  referenced.  Ths  earth  plane  is  designated  the  x,  y  plane  and  in  this  con¬ 
text,  the  antenna  coordinates  are  (0,  0,  c,).  The  aircraft  is  at  an  altitude 


z2  above  the  plane  and  at  position  and  on  the  plane.  Thus,  the  aircraft's 
position  is  (x2>  y2>  z^).  The  range  from  the  antenna  is  given  by  where 

RT=  f(z2  -  zi)2tx^y22]l/2  •  (2.1) 

This  geometry  is  shown  in  Fig.  1.  Another  important  constant  is 
the  ground  range  R  where 

R=[*2  +  y2J1/2  • 

Actually,  the  earth's  surface  is  not  a  flat  plane  but  contains  hills  and 
other  large  obstructions  such  as  buildings  and  bridges.  These  may  be  con¬ 
sidered  large  scale  surface  perturbations,  being  the  size  of  many  wavelengths 
of  the  transmitted  radiation.  There  are  also  small  scale  obstructions  such 
as  3mall  rocks,  trees,  grass,  water,  roadway  surfaces,  etc.  which  are  in 
size  less  than  several  wavelengths.  It  is  necessary  to  consider  both  of 
these  classes  of  obstructions  in  any  model.  In  many  cases,  the  large  scale 
obstructions  are  limited  in  number  and  are  easily  defined  and  characterized. 
For  example,  a  large  hill  can  be  easily  identified  from  a  topographic  map. 
However,  the  small  scale  obstructions  are  quite  numerous  and  detailed  in 
their  shape  as  compared  to  a  long  smooth  hill.  Thus,  the  large  scale  obstruc¬ 
tions  are  often  amenable  to  a  deterministic  analysis  whereas  random  analyses 
are  necessary  for  the  understanding  of  the  small  scale  effects.  Thus,  we 
shall  leave  discussion  of  the  latter  to  Section  4  and  consider  the  former  here. 


Fig.  1.  Geometry  of  canonical  coordinate  system 
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We  shall  assume  that  any  large  scale  obstructions  can  be  represented 
as  a  suitably  sized  section  of  a  plane  that  has  been  positioned  to  coincide  with 
the  actual  obstruction.  For  example,  if  there  is  a  large  building  parallel  to 
the  y  axis  in  Fig.  1,  then  it  can  be  represented  by  a  plane  of  equivalent  size 
and  electrical  properties  located  at  the  same  position.  These  equivalent 
planes  will  be  defined  relative  to  the  canonical  coordinate  system. 

A  model  of  the  scattering  surface  is  constructed  in  the  fashion  shown 
in  Fig.  2.  From  the  canonical  system  we  displace  another  orthogonal  coordi¬ 
nate  system  (x',  y1,  z')  by  a  vector  £q  (Fig.  2(a)).  Then  this  system  is  rotated 
about  the  y'  axis  by  an  angle  a  (see  Fig.  2(b))  called  the  cross -range  tilt.  It 
is  called  cross -range  because  the  y-axis  is  usually  associated  with  the  hori¬ 
zontal  range  and  multipath  reflections  from  such  a  tilted  surface  will  cause 
monopulse  azimuth  errors  in  the  cross -range  direction.  The  resulting  coor¬ 
dinate  system  x",  y",  z"  (Fig.  2(b))  is  then  rotated  about  x"  to  an  angle  0 
(see  Fig.  2(c))  yielding  a  final  coordinate  system  x'",  y1",  z'".  Miscalled 
the  range  tilt  angle.  Finally,  in  this  coordinate  system  on  the  x"',  y'"  plane, 
we  define  a  surface  S  (see  Fig.  2(d))  as  the  scattering  surface.  In  this  fashion 
we  can  construct  a  set  of  contiguous  surfaces  which  represent  the  large  scale 
surface  topography  and  present  a  one-to-one  mapping  of  the  large  scale  sur¬ 
face  topography.  In  the  next  section  we  shall  particularize  this  further  in  the 
study  of  specular  reflections.  We  shall  also  see  in  Section  4  that  surface 
roughness,  namely  small  scale  surface  irregularities,  can  be  included 
directly  by  superposition. 


(b)  Cross -Range  Tilt 


Fig.  Z.  Scattering  plane  goemetrv, 
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(C)  Ronge  Tilt 
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<j)  Scattering  Aperture 


tt taring  plane  goemetry  (Continued) 
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2.2 


SIGNAL  MODEL 


Now,  with  the  description  of  the  geometry  of  the  environment,  it  is 
possible  to  discuss  the  nature  of  the  signal  propagation.  We  shall  concen¬ 
trate  on  the  signal  transmitted  by  the  aircraft  and  received  by  the  sensor 
antenna,  which  we  term  the  downlink  signal.  For  large  scale  irregularities, 
such  as  hills,  we  cau  use  the  simplifying  analytical  techniques  of  ray  theory 
to  determine  the  propagation  of  the  electromagnetic  energy  from  the  transmit¬ 
ter  to  receiver.  For  example,  in  Fig.  3,  we  depict  the  aircraft  at  T  and  the 
receiver  at  R.  We  assume  that  the  x,  y  plane  is  the  surface  of  the  flat  earth 
model  and  that  S  is  -  separate  large  scale  scattering  surface.  The  points  P 
and  G  are  the  reflection  points  on  these  two  surfaces.  Path  (T)  is  called  the 
direct  propagation  path  and  represents  the  free  space  propagation  from  T  to 
R.  Path  (2)  ,  TGR,  is  the  ground  propagation  path  with  reflection  point  G. 

The  reflection  point  is  determined  directly  from  the  image  antenna  at  R', 

Path  (3),  TPR,  is  what  we  shall  call  a  multipath  signal  return.  It  is  a  sig¬ 
nal  return  that  results  from  other  than  flat  earth  reflections.  It  is  this  type 
of  return  that  proves  to  be  most  detrimental  to  direction  finding  (DF)  accuracy. 

The  signal  returns  in  Fig.  3,  are  clear,  well  defined  signals  which 
appear  at  the  receiver  as  signals  coming  from  directions  other  than  the  direct 
propagation  path.  Such  strong  well  defined  returns  are  called  specular  mul¬ 
tipath  and  result  from  large  reflecting  surfaces  such  as  buildings  and  hills. 

A  second  kind  of  multipath  signal,  called  diffuse  multipath,  occurs  due  to  the 
small  scale  surface  irregularities.  Obstructions  such  as  trees,  windows, 
rocks,  and  other  "small"  obstructions  cause  multipath  returns  to  appear  to 
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come  from  a  spectrum  of  directions  instead  of  a  single  well  defined  direction. 
An  example  of  such  a  phenomenon  is  shown  in  Fig.  4  where  we  depict  many 
small  reflectors.  Diffuse  multipath  is  a  random  phenomenon  as  a  result  of 
large  number  of  objects  giving  rise  to  it.  As  such,  an  analysis  of  its  effect 
must  be  phrased  in  terms  of  these  random  orientations.  Furthermore,  the 
analysis  of  diffuse  multipath  effects  requires  the  inclusion  of  diffraction 
phenomena  and  cannot  be  performed  using  classical  ray  optics  techniques. 

To  observe  the  effects  that  these  different  phenomena  have  on  the 
received  signal,  a  model  for  such  a  signal  must  be  developed.  In  general, 
the  signal  will  be  narrowband,  centered  about  a  carrier  frequency  uc>  We 
assume  that  the  transmitter  emits  a  spherical  wave  from  an  omnidirectional 
antenna  (variations  from  which  will  be  considered  later)  having  a  time  varia¬ 
tion  of  the  form. 


s(t)  =  J2.  Re[s(t)  exp(jwct)]  . 


(2.3) 


The  tevm  s(t)  is  the  complex  envelope  of  the  signal  and  has  the  form 


s(t)  *  ^  f(t) 


(2.4) 


where  E  is  the  signal  energy  and  f  (t)  is  a  normalized  version  of  the  temporal 
s 

variation.  Namely, 


f(t)  f  *(t)  dt  =  1 


(2.5) 
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Fig.  4.  Diffuse  multipath  reflection. 


The  signal  f  (t)  is  the  complex  envelope  of  the  time  variation  and  can 
take  on  any  arbitrary  form  depending  on  the  signal  format  required,  The 
signal  f  (t)  corresponds  to  the  video  output  (e.g.,  baseband  signal,  see  Van 
Trees). 

We  now  want  to  consider  the  complex  envelope  of  the  RF  signal  at 
the  output  of  the  antenna.  This  signal  will  be  affected  by  the  transmitter, 
the  transmission  channel,  and  the  antenna.  To  understand  the  antenna  effects 
we  must  first  define  the  incident  characteristics  of  a  plane  wave,  via  its  ray, 
upon  the  antenna. 

/  In  Fig.  5  we  have  drawn  an  array  in  the  canonical  coordinate  system 

and  have  defined  a  normal  vector  n  relative  to  the  aperture.  An  incident  plane 
wave  is  described  by  a  wave  vector  k.  The  elevation  angle  of  arrival  of 
the  plane  wave,  <J»,  is  the  angle  between  the  wave  vector  and  the  plane  formed 
by  the  vector  n  and  the  x  axis.  The  azimuth  0  is  defined  by  the  angle  made  in 
the  n,  x  plane  as  shown.  The  antenna  attenuates  the  amplitude  of  a  signal 
coming  from  direction  4>.  0  by  a  factor  G(0,  <|>)  which  is  called  the  antenna 
gain  in  that  direction.  Thus,  the  envelope  of  the  RF  signal  for  one  particular 
path  can  be  given  by 

G(0,  4>)  £  f  (t  -  R/c)  (2. 6) 

where  R  is  the  path  length,  ip  the  path  phase,  A  the  area  of  the  antenna,  p  the 
channel  transmission  factor  and  c  is  the  propagation  velocity. 

A  general  model  for  the  RF  envelope  can  now  be  proposed.  If  there 
is  a  free  space  signal  and  N  specular  multipath  signals  incident  on  the  antenna, 
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Fig.  5.  Definition  of  azimuth  and  elevation  angles. 
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plus  diffuse  multipath  and  noise,  then  the  received  signal  is  given  by; 

N 

r{t)  =  jEaA'y^j  G(QV  (jj^)  f  (t  -  Rj/ c)  exp(j^) 
i=0  1 


G(0,  <j>)  b(0,  cj)yt)  d0d<f>  +  n(t) 


(2.7) 


Here  0^  ^  are  the  azimuth  and  elevation  of  the  ith  specular  signal  and  Ri 
and  ^  are  the  transmission  gain,  range  and  phase  respectively.  The  i=0  terms 
correspond  to  the  free  space  path.  The  term  b(0,  <j>,  t)  represents  the  effects 
of  diffuse  multipath  to  be  discussed  in  Section  4.  The  noise  n(t)  represents 
all  noise  not  attributable  to  signal  transmission  observed  at  the  RF  output. 

One  of  the  purposes  of  this  report  is  to  describe  how  one  can  evaluate 
all  of  the  terms  that  appear  in  this  expression  and  show  how  they  relate  to 
the  actual  terrain  topography.  For  example,  the  quantity  p  will  depend  on 
reflection  coefficients,  diffraction  effects,  diffuse  reflection  effects,  as  well 
as  the  vagaries  of  the  transmitting  antenna.  The  term  b(0,  <$>,  t)  is  a  random 
process  imbedded  in  a  random  field.  The  evaluation  of  its  statistics  is  quite 
involved  and  will  require  some  detailed  knowledge  of  the  surface  topography. 


SECTION  THREE 


SPECULAR  MULTIPATH 

The  dominant  multipath  effect  in  terms  of  direction  finding  and  com¬ 
munications  in  the  DABS  environment  is  specular  multipath.  As  we  shall  see 
in  the  next  section,  diffuse  multip  ith  plays  a  second  order  role  which,  though 
marginally  important,  does  not  limit  the  performance  of  the  DABS  system. 

Specular  multipath  is  felt  in  terms  of  signal  fading  and  azimuth  estimation 
errors.  It  can  be  analyzed  deterministically  for  each  specific  geometry  or 
probabilistically  for  an  ensemble  of  such  geometries.  It  has  been  found  that 

i 

a  deterministic  analysis  provides  greater  insight  into  how  it  can  affect  system 
performance.  In  this  section  we  shall  concentrate  on  a  deterministic  analysis 
and  briefly  comment  on  proposed  random  models. 

As  introduced  in  the  last  section,  the  direct  xree  space  signal  pluse 
the  specular  signal  received  at  the  RF  output  of  the  antenna  is  given  by 

Jea  £c<V  *i>  IC  '  (* '  -r)  “POM  (3- l) 

i=0 

where  i=0  corresponds  to  the  direct  path  ray.  The  problem  is  to  determinia-  *  J 

| 

tically  evaluate  8^  4^,  R^,  and  ^  for  each  of  the  rays.  This  is  to  be  done 

knowing  the  number  of  scatterera,  the  surface  geometry,  the  antenna  pattern,  * 

the  aircraft  position  and  the  electrical  properties  of  the  surface. 
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3.1 


GEOMETRICAL.  CONSIDERATIONS 


To  evaluate  these  quantities,  we  must  first  further  quantify  the  geo¬ 
metry.  Let  us  assume  that  the  receiving  antenna  is  an  array  and  the  phase 
(geometrical)  center  of  the  array  is  at  some  vector  r  ^  where 


Il  = 


X1 

*1 

Z1 


(3.2) 


which  is  a  generalization  of  the  definition  of  Section  Two.  Likewise,  assume 
the  aircraft  is  at  r ^  where 


?2 


(3.3) 


l2 


2J 


Both  of  these  vectors  are  defined  relative  to  the  canonical  coordinate  system 
developed  in  the  last  section.  Now  we  will  <  .  "'ve  the  ith  scattering  plane  (i  £  2) 
as  the  septuple  of  surface  area  Sj,  with  tilts  |3.(i  £  2)  where; 


(3.4) 


Here,  r^  is  the  offset  vector  of  the  ith  scattering  plane  (see  Fig.  2),  ^  and 

t r .  are  the  dielectric  constant  and  conductivity  if  the  plane.  <r  is  the  surface 
i  si 

roughness  coefficient  as  obtained  for  diffuse  multipath  effects.  This  septuple 


completely  defines  the  scattering  plane  and  will  provide  all  the  information 


necessary  for  our  model  of  the  channel.  Given  the  antenna  gain,  the  received 
signal  is  completely  defined. 

We  shall  proceed  by  first  determining  R.,  the  range  between  the  trans¬ 
mitter  on  the  aircraft  and  the  receiver  via  the  ith  multipath  source.  We  shall 
briefly  outline  this  derivation  and  present  the  result.  When  using  a  ray  optics 
approach  the  point  of  reflection  from  the  ith  surface  is  determined  by  extend¬ 
ing  the  ray  from  the  source  to  the  image  antenna.  The  image  antenna  is  de¬ 
fined  as  being  that  point  relative  to  the  scattering  plane  (xm ,  y"1  plane)  which 
is  equidistant  from  the  opposite  side  of  the  plane  as  the  receiver  and  lying  on 
the  normal  created  by  the  receiver  and  the  scattering  plane.  For  example,  in 
Fig.  6  we  have  depicted  an  antenna  at  R  in  the  y,  z  plane  and  a  transmitter  at 
T  also  in  the  y,  z  plane.  The  image  antenna,  since  the  x,  y  plane  is  the  scat¬ 
tering  pi  me,  is  at  the  point  R.  The  line  RO  is  normal  to  the  x,  y  plane. 

The  length  of  the  path  is  T5  +  SR  or  equivalently  TR1.  The  vector  r^  is  the 
vector  from  the  origin  to  the  image  antenna.  What  should  be  noted  in  this 
figure  is  that  since  the  x,  y  plane  is  the  scattering  plane  and  the  image 
antenna  lies  along  the  normal  RO  at  the  same  distance  from  the  x,  y  plane 
as  rj,  then  has  the  same  x,  y  coordinates  are  £j  and  a  z  coordinate 
that  is  the  negative  of  rj.  This  then  suggests  how  the  image  antenna  for  the 
ith  scattering  plane  can  be  obtained.  Namely,  translate  and  rotate  the  coor¬ 
dinate  system,  relative  to  the  canonical  system,  using  r^  ,  and  and  em¬ 
ploy  the  algorithm  suggested  above.  Then  retranslate  and  rotate  back  to  the 
canonical  coordinate  system.  Doing  this  we  find  that  the  position  of  the  image 
antenna  in  the  canonical  coordinate  system  is  given  by  r',  where 


(3.5) 


£il='2i<£i  -£o,)  +  I- 


1 


where  is  the  matrix  which  performs  the  above  mentioned  operations.  The 
entries  of  &  are; 


1 


Si 


2  2 

cos  -  sin  a.  cos  2^ 

since  sin2$. 

2 

sin2a*cos  0. 

sina.  sin20 

cos20. 

cosa.sin2|3. 

i  i 

l 

x  l 

-sin2a^  cos^3^ 

cosaj  sin2/3^ 

2  2 

sin  -  cos  a.cos2P| 

(3.  i ) 


Thus  R^,  the  path  length,  via  the  ith  multipath  reflector,  is  deter¬ 


mined  by 


1  )T(r,  -  r\  H1/2 


Ri  ■  [<£z  '  £)t>  <£z  -  li.'l 


(3.7) 


This  is  the  value  of  R^  used  in  (2.7).  Now  ^  represents  the  phase  of 
the  ith  path.  If  we  arbitrarily  define  then 


*i**0 


+  2ir 


(Rt- 


V 


F. 

l 


(3.8) 


where  is  the  free  space  wavelength  and  Rq  is 


The  term  ^  is  the  phase  of  the  ith  reflection  process  which  we  shall  define 
i 

shortly  (see  3.  2.  1). 

In  a  similar  fashion,  we  can  use  this  formalism  to  evaluate  the  azi¬ 
muth  and  elevation  (9^,  4>^)  of  the  ith  multipath  source.  Consider  the  geome¬ 
try  in  Fig.  7(a).  The  xy  plane  is  the  reflection  plane  which  contains  a  point  of 
reflection,  S.  Note  that  at  the  point  S  the  ray  changes  direction  and  is  reflected 
to  the  antenna.  The  azimuth  and  elevation  are  then  defined  relative  to  the 

vector  that  represents  the  scattered  wave.  Let  us  call  that  vector  r^  .  Now 

i 

r  relative  to  an  arbitrary  shift  of  coordinates  is  equal  to  the  vector  the 
— 8i 

incident  vector,  except  that  the  z  coordinate  has  its  sign  changed.  From  the 
geometry,  is  given  by 


£i  =  *2  ■  £it  •  (3- 10) 


This  is  defined  as  the  vector 


(3.11) 


Consider  Fig.  7(b).  Here  we  depict  r  ,  the  scattered  version  of  s..  The 
elevation  of  the  multipath  relative  to  the  array  is  defined  as  -  4>j,  as  shown, 
and  the  azimuth  is  -  0^  as  shown.  Thus  and  0^  are  easily  obtained  by  know¬ 
ing  the  components  of  r  which  are  directly  obtained  from  s.  as  described. 

This  technique  can  then  be  extended  to  any  arbitrarily  rotated  plane 

following  the  reasoning  developed  before.  Namely,  r  can  be  shown  to  be; 

”®i 
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ay  i  '3iTT  ****  fWt 


(3.12) 


1 


so  that  the  azimuth  is  given  by 


=  tan 


-1 


*il 


i2j 


(3.13) 


and  the  elevation; 


4>i  =  tan 


-1 


s 


i3 


,  i — :~n 

IV  r‘ 


i2 


T/2 


where 


{ 

(3. 14) 


i 


(3.16) 

1 


Thus,  given  the  parameters,  we  can  obtain  8^,  <J>^  and  then  knowing  the  form  of 
the  antenna  gain  G(6j,  <$>^)  we  obtained  the  antenna  dependent  portion  of  the  ith 
return. 

Another  term  that  is  of  interest  is  the  angle  made  between  the  scatter¬ 
ing  plane  and  the  incident  ray.  As  we  shall  see,  this  is  important  in  evaluating 
the  Fresnel  reflection  coefficient  of  the  scattering  surface.  In  the  case  where 
the  scattering  plane  is  the  x,  y  plane  then  as  shown  in  Fig.  8  this  angle  is 
given  by  the  angle  ?.  The  angle  TSR,  equals  v  -  2£.  Now  consider  the 
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general  case.  Let  v 3  be  a  vector  from  the  origin  to  the  receiver  at  R  and  let 
Vj  be  a  vector  to  the  image  receiver  at  R'.  Then  the  vector  v,  -  v^  =  v^  is 
normal  to  the  plane  of  scattering.  For  example,  in  Fig.  9  it  is  normal  to  the 
x,  y  plane.  Now  let  v^  be  a  vector  to  the  transmitter.  Let  Vg  be  the  vector 
from  the  transmitter  to  image  antenna.  Now  this  vector  is  in  the  direction 
of  the  ray  incident  on  the  scattering  plane  for  which  v^  is  a  normal  vector. 
Let  y  be  the  angle  between  v^  and  Vg.  This  angle  can  be  obtained  from  the 
relationship. 


Then  by  this  simple  geometrical  argument  ?  equals  ((ir/2)  -  y)  (see  Fig*  9). 

3.  2  REFLECTION  COEFFICIENT 

The  last  factor  required  for  (3. 1)  is  an  evaluation  of  the  reflection 
coefficient,  p^.  This  coefficient  depends  upon  the  electrical  properties  of  the 
surface,  the  size  of  the  surface,  and  the  small  scale  surface  perturbations. 
The  electrical  properties  are  used  to  evaluate  the  Fresnel  reflection  coeffi¬ 
cient,  which  has  both  an  amplitude  and  phase.  The  size  of  the  scattering  sur¬ 
face  defines  an  aperture  and  the  surface  roughness  determines  the  diffractive 
effects  that  may  arise.  This  is  a  very  complex  issue  and  must  be  dealt  with 
in  detail.  Finally,  diffuse  scattering  behavior  of  the  large  scale  surface  may 
tend  to  disperse  the  incident  energy. 
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Considering  the  above  three  effects  the  reflection  coefficient,  p^, 
can  be  written  as; 


pi  =  51  A{  r1Di  rlFi 


(3. 17) 


where  r]A  is  the  coefficient  accounting  for  the  finite  aperture  effect,  rjn  the 

1  u{ 

diffuse  reflection  coefficient  and  r the  Fresnel  reflection  coefficient  due  to 

Fi 

the  electrical  properties  of  the  surface. 


3. 2. 1  Fresnel  Reflection  Coefficient 


The  Fresnel  reflection  coefficient  rjp  and  the  phase  <p-p  ,  (3.8),  can 
be  obtained  by  an  analysis  of  the  reflection  of  electromagnetic  waves  from 
smooth  surfaces  with  known  electrical  properties.  The  properties  of  the 
reflector  also  depend  upon  the  polarization  of  the  incident  plane  wave.  To 
simplify  the  analysis  we  shall  assume  a  single  polarization,  vertical,  and 
avoid  until  Section  5  any  of  the  problems  associated  with  depolarizing  surfaces. 

The  complete  reflection  coefficient  is  a  complex  quantity  with 


r,Fl  =  r,Fl  exP^^Fi) 


(3.18) 


It  is  shown  In  Jordan  and  Balmain  (p.  631)  that  for  vertical  polarization 


^  (er  -  jx)  -  \J(et  -  jx)  -  cos2? 
*i  (er  -  jx)  sin?  +  \/(er  -  jx)  -  coa£? 


(3.19) 
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where  j  equals  -J- 1  and  £  is  the  angle  of  incidence  given  by  (ir/2)  -  y,  where 
Y  is  defined  in  (3. 16).  x  is 


X  =  <7/U€v 


(3.20) 


where  <r  is  the  reflector's  conductivity,  u  the  angular  frequency  of  the  inci¬ 
dent  radiation  and  the  dielectric  constant  of  the  propagation  medium.  Also* 


er  =  e/e, 


(3.21) 


where  e  is  the  dielectric  constant  of  the  scattering  surface. 

Thus  for  different  electrical  surface  constants,  e  and  £ r,  values  of 

and  can  be  obtained  as  a  function  of  the  incident  angle  £.  Usually  £ 

Fi  Fi 

is  quite  small,  which  means  that  rj_  is  close  to  unity  and  close  to  it. 

*i  i 

However,  for  different  surfaces,  the  behavior  of  these  terms  as  a  function  of 
£  varies  significantly.  In  Fig.  10  we  have  shown  the  magnitude  and  phase  for 
a  typical  surface.  It  can  be  seen  from  Fig.  10  that  a  significant  change  in 
r/p  occurs  as  £  Increases.  Also  the  dip  at  the  Brewster  angle  may  vary  from 
a  significant  drop  to  only  a  3  IB  drop  in  amplitude  depending  upon  the  surface 
characteristics . 


3.  2. 2  Diffuse  Reflection  Coefficient 

When  the  surface  of  the  scatterer  is  rough,  energy  from  the  incident 
ray  is  scattered  away  from  the  specular  scattering  direction  (and  thus  is  lost 
from  the  received  signal).  To  account  for  this  effect,  we  introduce  the  diffuse 
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Fig.  10.  Reflection  coefficient 
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reflection  coefficient  r)n  *or  each  surface.  This  coefficient  has  been  deter - 
mined  by  Beckmann  [2]  and  it  depends  upon  the  angle  of  reflection,  £,  as  does 
the  Fresnel  reflection  coefficient.  It  also  depends  upon  the  statistical  proper¬ 
ties  of  the  scattering  surface.  Specifically,  if  z(x,  y),  the  height  of  the  sur¬ 
face  at  point  x,  y,  is  a  zero  mean  Gaussian  random  variable  with  variance 
2 

o^.  Beckmann  shows  that 


/4ir<rh  sin  A  2 

n°i=  eXP[  \  ^  /  ‘ 


(3.  22) 


For  ?  close  to  0  the  coefficient  is  near  unity.  However,  for  very  rough 
surfaces  (tr^  >>  \q),  this  term  decreases  rapidly  as  the  angle  of  inci¬ 
dence  increases.  In  Fig.  11  we  have  sketched  r}n  for  two  deferent  roughness 

l 

ratios,  °^Aq.  as  a  function  of  5. 

3. 2.  3  Diffraction  Effects 

As  the  size  of  the  reflecting  surface  decreases,  diffraction  becomes 
important  a.nd  ultimately  dominates  the  behavior.  For  example,  it  is  well 
known  that  a  rectangular  surface  illuminated  by  a  plane  wave  from  direction 
k  as  shown  in  Fig.  12(a)  diffracts  the  radiation  into  the  specular  direction  and 
into  other  directions  according  to  the  sin  x/x  distribution.  In  Fig.  12(b)  we 
have  plotted  the  distributions  of  the  diffracted  amplitude  of  the  field  in  the  k^, 
ky  directions.  Note  the  central  peak  in  the  specular  direction,  k\  and  the 
presence  of  sidelobes  in  other  directions.  In  Fig.  13  two  representative 
specular  reflectors,  and  S2<  are  shown  with  sketches  of  the  diffraction 
pattern  of  the  reflected  radiation  superimposed.  The  dimensions  of  surfaces 


>  * 
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S  j  and  are  assumed  to  be  small  enough  that  the  radius  of  curvature  of  the 
wavefront  is  large  compared  to  the  surface  size.  This  allows  us  to  say  that 
the  incident  wave  is  coming  from  the  direction  joining  the  midpoint  of  the  sur¬ 
face  to  T,  and  the  the  surface  is  irradiated  by  a  plane  wave. 

The  amplitude  of  the  field  received  from  a  reflector  depends  on  the 
width  and  orientation  of  the  diffraction  pattern,  which  in  turn  depends  on  the 
size  and  location/orientation  of  the  reflector,  respectively. 

From  geometrical  optics  we  know  that  if  T  emits  a  spherical  wave 
and  the  wave  is  scattered  by  some  plane,  then  the  scattered  wave  is  also 
spherical.  Furthermore,  the  apparent  direction  o£  the  radiation  can  be 
obtained  by  observing  the  direction  of  the  ray  passing  through  the  desired 
point.  Thus,  in  Fig.  14,  we  have  drawn  a  ray  diagram  for  a  wave  travelling 
from  T  and  have  shown  seven  receivers.  The  reflected  wave  is  also  shown  to 
be  spherical  by  drawing  the  normal  to  the  rays.  Now  consider  any  one  of  the 
receiver  points  R^,  For  this  case  of  an  inifinite  reflecting  plane,  we  see  that 
there  is  a  definite  reflecting  point,  P.,  at  which  the  ray  follows  Snells  law. 

We  now  pose  the  question:  If  there  is  only  a  finite  amount  of  scattering  sur¬ 
face,  and  Pj  is  not  on  that  surface,  what  is  received  at  R^?  To  find  the 
answer  we  must  consider  diffraction  analysis.  Thus,  in  Fig.  13,  if  R  is  in 
the  sidelobes  of  that  diffraction  pattern  of  surface  Sj,  its  amplitude  is  de¬ 
creased.  Furthermore,  we  can  consider  the  following  experiment.  If  we 
fix  R  and  T  in  Fig.  13,  move  the  surface  S  and  plot  the  amplitude  of  the 
received  field  as  a  function  of  y,  the  range,  we  will  effectively  plot  out  the 
diffraction  pattern  for  that  surface.  This  result  is  sketched  in  Fig.  15.  When 
the  surface  is  centered  at  a  point  y*,  the  amplitude  of  the  signal  received  at 


R  reaches  a  maximum 


GROUND  PLANE 


Now  if  S  is  large  enough  (to  have  a  narrow  diffraction  pattern)  but 
not  too  large  (small  with  respect  to  the  radius  of  curvature)  then  the  point  y* 
in  Fig.  15  will,  in  some  sense,  correspond  to  the  reflection  points  Pj  in 
Fig.  14.  That  is,  by  positioning  this  small  surface  at  the  right  point,  we 
can  effect  a  received  wave  that  would  appear  to  be  coming  from  an  infinite 
reflecting  plane.  This  idea  is  the  basis  of  the  Fresnel  zone  concept.  Further 
more,  it  can  be  rigorously  shown  that  if  one  has  many  such  planes,  that  only 
the  one  centered  at  the  right  position  will  contribute  to  the  received  field 
while  all  others  will  coherently  cancel. 

In  view  of  the  above  discussion,  we  can  consider  what  portion  of  an 
infinite  reflecting  plane  dominates  in  the  reflecting  process.  This  area  is 
called  the  first  Fresnel  zone  and  represents  the  region  on  the  plane  where 
the  first  ir  radians  of  phase  shift  occurs.  The  higher  order  Fresnel  zones 
represent  regions  of  increasing  rr  radians  of  phase  shift.  These  zones  are 
depicted  in  Fig.  16(a).  The  first  zone  is  shown  in  Fig.  16(b)  in  more  detail. 

It  is  an  elliptical  surface  located  at  a  center  point  y^j  with  an  extent  +  A  j 
about  that  point  and  +  Xj  about  the  x  axis.  These  quantities  are  given  in 
terms  of  Zj,  z^,  and  r  as  follows: 


*01  8 


1  + 


Zzjfzj  +  z2yi 


1  + 


(z.  +  z?) 


r 

2 


(3.23) 


/ 

39 


26J1/2 


a  -r  L  "OJ 
Ai "  z  — ; — t — 7 

i  +  -A-r- 

V 


(3.24) 


where  6q  is  defined  by; 
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The  area  of  the  Fresnel  zone,  A^,  equals; 


af=  it  XjAj  . 


(3.27) 


With  these  values  we  can  determine  whether  or  not  a  surface  will  act  as  a 
source  of  reflection  from  the  transmitter  to  receiver.  Namely,  if  it  is  in  the 
first  Fresnel  zone  then  it  clearly  will  reflect  an  amount  limited  by  the  dilution 
of  the  main  lobe  of  the  diffraction  pattern  due  to  its  finite  size,  If,  however, 
it  lies  ouside  the  first  Fresnel  zone  then  its  contribution  to  the  reflected  sig¬ 
nal  will  be  limited  to  that  portion  which  is  in  the  zone.  Thus,  as  the  reflect¬ 
ing  surface  moves  out  of  the  Fresnel  zone  its  reflected  signal  amplitude  at 
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the  receiver  decreases,  the  amplitude  being  determined  by  the  sidelobes  of 
the  diffraction  pattern. 


Thus  an  evaluation  of  the  diffraction  reflection  coefficient,  rjn  ,  can 

ui 

be  determined  from  the  diffraction  pattern  of  the  surface  being  irradiated. 
However,  such  a  calculation  is  often  quite  tedious.  A  reasonable  approxima¬ 
tion  is  let  rjn  equal  the  ratio  of  the  amount  of  reflecting  surface  area  in  the 
first  Fresnel  zone  to  the  area  of  the  first  a'resnel  zone.  This  is  a  quantity 
that  is  easily  calculated.  That  is,  if  is  the  area  of  the  ith  reflector  in  the 
first  Fresnel  zone  then; 
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Thus,  (3.19),  (3. 22)  and  (3.  28'  provide  all  that  is  necessary  for 
in  (3.  17).  Furthermore,  this  completes  the  specifications  of  all  the  para¬ 
meters  in  the  specular  signal  of  (3.  1). 

3.  3  EXAMPLES  OF  RECEIVED  SIGNALS 

The  model  developed  at  the  beginning  of  this  section  for  the  specular 
portion  of  the  received  signal  can  be  used  to  demonstrate  the  effect  of  certain 
terrains  or  structures  on  the  signal.  The  dominant  effect  is  that  of  fading, 
where  the  direct  return  and  other  multipath  returns  add  coherently  to  cause 
a  decrease  in  signal  amplitude.  This  can  produce  a  reduction  in  the  signal  - 
to-noise  ratio  which  seriously  affects  detection  and  position  estimation.  Using 
the  model  just  developed,  it  is  possible  to  see  what  type  o*  multipath  environ¬ 
ment  will  give  rise  to  the  more  deleterious  effects  and  how  through  proper 
siting  and  antenna  design,  they  may  be  prevented. 
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The  fading  property  of  the  signal  is  easily  observed  by  evaluating  the 
power  in  the  specular  returns.  Thus,  let  <r(t)  be  the  magnitude  squared  of 
(3. 1).  It  can  be  written  as; 
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Now  if  f(t)  has  a  width  which  is  long  in  time  with  respect  to  the  path  differences 
then  the  time  dependence  may  be  neglected  by  integrating  (3.  29)  and  using  the 
assumption  in  (2.  5).  The  resulting  time  integrated  function  is  called  s  and  is 
(3.  29)  with  the  time  dependence  absent.  The  first  terms  are  clearly  indepen¬ 
dent  of  the  path  differences  and  are  always  positive.  The  second  terms,  those 
involving  the  double  summation,  depend  on  the  path  differences  due  to  ^  -  <//. 
and  can  be  negative.  They  represent  the  result  of  coherent  addition  of  the 
wavefronts. 

This  phenomenon  has  been  investigated  by  Spingler  and  Fig.  17  is  a 
plot  of  these  amplitudes  versus  aircraft  elevation  angle  on  a  radial  flight. 

One  can  note  clearly  the  beat  phenomenon  and  see  the  multiple  beating.  These 
results  are  for  a  typical  ATC  radar  beacon  interrogator  antenna.  Fades  of 
20  dB  are  not  uncommon  in  this  data  set. 
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Using  the  model  developed  in  this  section  we  have  analyzed  three 
specific  cases  in  detail.  Each  case  assumed  a  flat  earth  plus  a  single  addi¬ 
tional  scattering  plane.  The  three  cases  are  depicted  in  Fig.  18.  Case  I  is 
an  example  of  a  slightly  sloping  hill.  Case  II  is  a  large  vetical  obstruction 
that  is  parallel  to  the  y  axis.  Case  HI  represents  a  moderately  sloping  rise. 
In  all  cases,  the  aircraft  flew  a  radial  out  along  the  y  axis  at  a  fixed  altitude. 

With  each  of  these  examples  we  also  analyzed  the  effect  of  multipath 
on  monopulse  azimuth  estimation  and  the  performance  of  multipath -inter¬ 
ference  detectors  (MID)  (see  McAulay  and  McGarty).  For  monopulse  we 
have  two  responses  2  and  A  corresponding  to  an  antenna  gain  Gs(8,  cj>)  or 
G^(0,  <)>)  respectively.  Using  these  gains  in  (3.  1)  and  properly  eliminating 
the  time  behavior  (matched  filter  detection)  one  obtains  2  and  A.  It  is  clear 
that  they  are  complex  numbers.  The  estimated  angle  is  given  by 


0  =  k  Re 


(3,  30) 


where  k  is  the  monopulse  slope.  The  MID  is  given  by  t  where 
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In  Figs.  19-21  we  have  plotted  s,  0  and  t  respectively  versus  l/R  (where  R 
is  range).  The  plot  of  s  in  this  fashion  shows  the  periodic  nature  of  the 
received  signal  directly.  The  rapid  beats  are  due  to  the  extra  reflecting 
plane  while  the  lower  frequency  modulation  is  due  to  the  flat  earth  reflection. 
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Fig.  19.  Beam  power  vs  inverse  range;  Case  I. 
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Fig.  ?.0.  Estimated  angle  vs  range;  Case  I. 
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Fig.  21.  Multipath  test  vs  range;  Case  I, 


The  angle  estimate  in  Fig.  20  for  Case  I  shows  the  combination  of  both  fading, 
SNR  losses,  and  the  multipath  bias  modulation  effect  note  by  McAulay  and 
McGarty.  The  plot  of  r  shows  how  such  multipath  is  registered  by  this  detec¬ 
tor  by  its  passage  beyond  the  detection  thresholds  shown.  These  cases  were 
done  for  512  element  dipole  array  that  was  tilted  down  4°  and  had  a  3  dB/degree 
vertical  cutoff  below  the  "horizon." 

In  a  similar  fashion,  Case  II  is  shown  in  Figs.  22-24.  Here  the  mul¬ 
tipath  effect  on  angle  estimates  is  more  clearly  pronounced  in  Fig.  23.  Case 
HI  results  are  shown  in  Figs.  25-27  for  a  hogtrough  antenna  and  in  Figs.  28- 
30  for  an  antenna  with  vertical  aperture.  What  is  most  clear  here  is  that 
vertical  aperture  does  help  significantly  in  reducing  the  angle  errors  (Fig.  29). 
This  is  a  result  of  increasing  the  effective  SNR  as  is  evidenced  by  comparing 
Figs.  2 5  and  28. 

Various  other  analyses  are  possible  that  exhibit  the  effects  of  in¬ 
creased  surface  roughness,  changes  in  the  electrical  properties  and  other 
reflecting  plane  orientations. 

3. 4  CONCLUSIONS 

The  effects  of  reflections  from  large  smooth  surfaces  has  been 
modeled  as  the  coherent  sum  of  plane  waves  incident  on  the  antenna  aperture. 
The  waves  can  be  completely  characterized  by  their  ampli  de,  phase, 
azimuth  and  elevation  relative  to  the  boresight  of  the  antenna.  In  ihis  section 
we  presented  a  detailed  model  which  could  be  used  to  develop  a  received  sig¬ 
nal  which  would  be  representative  of  specular  multipath. 

The  specular  multipath  model  developed  depended  upon  the  following 

points : 
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Fig.  27.  Multipath  detector  statistic;  Case  IH,  no  vertical  aperture 
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1. 


Geometry:  For  each  multipath  source  a  plane  was  posi¬ 
tioned  at  the  source  location  and  tilted  in  both  range  and 
cross  range  directions  to  coincide  with  the  actual  surface. 
Evaluation  of  propagation  distances  could  then  be  made 
directly  knowing  the  location  of  the  receiver  and  trans¬ 
mitter.  All  geometrical  constants  were  determined  for 
each  multipath  reflector  in  terms  of  the  offset  position  of 
the  plane,  its  range  and  cross  range  tilts  and  the  locations 
of  transmitter  and  receiver. 

Reflection  Coefficient:  The  effect  of  a  multipath  signal  on 
azimuth  estimation  and  link  reliability  depends  directly 
upon  how  strong  it  is  relative  to  the  free  space  signal. 

The  ratio  of  multipath  to  free  space  signal  amplitudes  is 
given  by  the  reflection  coefficient.  These  were 

a.  Electromagnetic  Reflection:  This  is  determined 
by  the  Fresnel  reflection  equation  for  a  plane 
wave  at  the  interface  of  two  media  with  differ¬ 
ent  electrical  properties. 

b.  Diffuse  Multipath:  This  represents  the  fraction  of 
incident  radiation  that  is  not  lost  to  other  than  the 


specular  direction  as  a  result  of  surface  roughness, 

c.  Diffraction:  If  the  surface  is  of  finite  size  than  its 


location  and  size  must  be  considered.  This  results 


in  a  study  of  where  the  scattering  aperture  falls 
relative  to  the  Fresnel  zone  projected  on  the  scatter 
ing  plane. 
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The  inclusion  of  all  these  items  leads  to  an  analytical  model  which 
can  be  used  to  simulate  what  would  be  observed  by  an  arbitrary  antenna  from 
a  given  terrain.  Specific  results  for  an  antenna  with  the  characteristics  of 
the  DABS  experimental  facility  antenna  (vertical  aperture  8X.q,  horizontal 
aperture,  26 XQ)  were  evaluated.  The  antenna  has  a  sum  beam  with  a  3.  2° 
beamwidth  in  azimuth  difference  beam  and  a  sector  beam  in  elevation  with 
3  dB/degree  cutoff  at  the  horizon.  Plots  of  the  sum  beam,  the  monopulse 
azimuth  estimate  and  the  multipath-interference  detector  were  plotted  for 
three  different  multipath  geometries.  As  a  result  of  these  simulations  the 
following  conclusions  were  reached. 

1.  The  monopulse  azimuth  errors  were  in  general  quite  small 
except  when  both  fading  from  the  flat  earth  and  the  multipath 
plane  (tilted  surface)  occured  simultaneously.  In  those 
instances  the  effective  free  space  signal  was  reduced  resulting 
in  an  exceptionally  large  multipath  to  signal  ratio  and  thus  a 
large  error.  This  effect  can  be  countered  however  by  using 

a  large  vertical  aperture  with  a  sharp  cutoff  below  the  hori¬ 
zon. 

2.  When  the  vertical  aperture  is  reduced  to  X^  serious  fades 
and  monopulse  errors  v, 

3.  The  inclusion  of  a  vertical  aperture  can  reduce  vertical  lobing 
fades  on  the  sum  beam  power  from  20  dB  with  a  Xq  aperture 
to  6-8  dB  with  an  8\^  aperture. 

4.  The  MID  shows  strong  correlation  between  large  azimuth 
errors  and  its  output.  Therefore,  it  may  be  useful  as  a  multi - 
path  detector. 
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The  simulations  discussed  in  this  report  represent  only  a  small 
number  of  those  that  have  been  performed.  Furthermore,  recent  compari¬ 
sons  of  simulations  with  actual  data  obtained  at  DABSEF  indicate  that  the 
model  does  represent  the  observed  phenonmena  quite  well  in  many  cases. 


SECTION  FOUR 


v 


DIFFUSE  MULTIPATH 


Diffuse  multipath  is  a  random  phenomenon,  distinct  from  specular 
multipath,  and  must  be  treated  as  such.  In  contrast  to  the  development  in 
the  last  section,  where  exact  descriptions  of  the  scattering  surface  were 
given,  in  this  section  the  nature  of  the  scattering  surface  is  describable  in 
only  a  probabilistic  sense.  The  reason  for  this  is  that  the  diffuse  multipath 
arises  from  the  roughness  of  the  surface  and  is  predominantly  a  diffraction 
effect.  This  will  become  important  when  we  analyze  the  effects  giving  rise 
to  the  scattered  field. 

In  the  estimation  of  the  azimuth  of  a  cooperative  signal  source,  the 
effects  of  multipath  tend  to  degrade  the  performance  of  such  estimates.  The 
multipath  encountered  has  been  divided  into  two  different  categories,  specular 
r  1  diffuse.  Specular  multipath  results  from  large  smooth  surfaces  occupy  ¬ 
ing  a  significant  part  of  the  first  Fresnel  zone.  Diffuse  multipath  diffeis 
from  the  specular  form  in  that  when  it  is  scattered,  it  does  not  propagate  in 
a  single  direction  but  in  a  continuum  of  directions,  depending  on  the  rough¬ 
ness  properties  of  the  scattering  surface.  Diffuse  multipath  is  a  diffraction 
phenomenon  where  the  element  doing  the  diffracting  is  a  small  surface  per- 
tubation. 

The  effects  of  such  scattering  were  discussed  by  Kerr(1951)  in  an 
attempt  to  analyze  the  effect  of  sea  surface  scattering  on  radar  performance. 
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The  approach  used  by  Kerr  was  to  obtain  a  scattering  cross  section  for  the 
surface.  A  more  detailed  approach  was  undertaken  by  Rice  (1931)  when  he 
modeled  the  surface  by  means  of  a  deterministic  polynomial  and  then  solved 
the  random  problem.  His  approach  was  qua  si -deterministic  and  directed  at 
sea  surface  scattering.  A  different  approach  using  the  Kirchoff  approxima¬ 
tion  was  presented  by  Eckart  (1953)  and  was  directed  at  acoustic  scattering 
from  rough  surfaces.  Equations  for  the  field  were  obtained  and  then  averaged. 
The  resulting  equations  were  then  solved  for  the  required  quantities.  This 
work  was  latter  followed  by  Ament  (1953,  1956),  using  a  similar  approach. 

Hoffman  (1955)  extended  the  analysis  to  the  vector  problem  encount¬ 
ered  in  electromagnetic  field  problems,  using  *he  Kirchoff  technique.  Other 
attempts  were  made  by  Clark  and  Hendry  (1964)  to  evaluate  backs cattering 
effects.  The  prime  reference  using  the  Kirchoff  method  is  Beckmann  and 
Spizzichino  (1963).  Their  results  give  the  field  decomposition  in  terms  of 
different  scattered  directions.  They  assume  a  Gaussian  surface  behavior 
with  knowledge  of  surface  standard  deviations  and  correlation  lengths.  Experi¬ 
ments  yielding  these  values  have  been  discussed  by  Fung  and  Moore  and  Hayre 
and  Moore. 

Most  of  these  previous  results  are  not  suitable  for  an  analysis  of 
azimuth  estimator  performance,  however.  It  was  shown  by  McGarty  (1974), 
for  the  case  of  bearing  estimation  with  a  sonar  array,  that  a  function  called 
the  channel  spread  function  was  necessary  to  evaluate  performance.  This 
function  gives  the  intensity  of  radiation  incident  at  the  receiver  from  all  direc¬ 
tions.  This  function  cannot  be  obtained  from  the  Beckmann  and  Spizzichino 
model.  However,  analyses  by  Kodis  (1966)  and  Barrick  (1968)  provide  exactly 
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what  is  necessary  for  the  performance  determination.  Their  analyses  obtain 

a  scattering  cross  section  which  relates  the  scattered  radiation  into  a  specific 

direction  from  a  specified  incident  direction.  The  method  of  analysis  differs  i 

greatly  from  others  in  that  it  solves  the  scattering  problem  first  using  a  ^ 

stationary  phase  technique  (similar  to  Twersky  (1957)),  and  then  an  averaging 

technique.  The  averaging  technique  yields  a  very  physical  formula  for  the 

scattering  cross-section.  It  is  in  terms  of  the  average  number  of  scatterers, 

the  average  surface  curvature  (following  Longuet -Higgins  [1,2])  and  reflection 

coefficient. 


4.  1  THE  SIGNAL  MODEL 

In  Section  2,  while  discussing  the  nature  of  the  total  returned  signal, 
we  defined  the  diffuse  multipath  return  as 


G(0,  4,)  b(0,  4> ,  t)  d0  dq> 


where  0  is  azimuth,  <j>  elevation  and  G(0,  <f>)  the  antenna  gain  for  those  angles. 
The  terms  b(0,  <j),  t)  represents  the  random  return  from,  0,  $  at  time  t.  It 
is  clear  then  that  b  (0,  <j>,  t)  is  some  form  of  a  random  process  in  time.  More¬ 
over,  it  is  also  a  random  process  in  space.  Such  processes  are  also  called 
random  fields. 

We  can  motivate  the  structure  of  b(0,  4>,  t)  by  considering  the  geome¬ 
try  shown  in  Fig.  31.  Here,  as  before,  we  have  a  transmitter  at  T  and  a 
receiver  at  R.  We  place,  centered  at  (x1,  y'),  an  incremental  surface  dS. 
Energy  is  now  transmitted  from  T  and  scattered  from  dS  into  all  directions. 
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The  amount  arriving  at  R  will  depend  on  the  diffraction  nature  of  the  surface  1 

(recall  our  discussion  in  Section  3.  2.  3).  \ 

i 

The  azimuth  0  and  elevation  <j>  are  shown  in  Fig.  31.  It  will  be  shown  i 

latter  that  each  (x1,  y')  pair  uniquely  corresponds  to  a  (0,  4>)  pair,  thus  with 
this  one-to-one  transformation  we  can  equivalently  phrase  the  multipath  in  i 

terms  of  these  variables.  Let  us  now  also  consider  the  temporal  behavior 
of  b(0,  <(>,  t).  Since  we  are  assuming  that  the  surface  is  not  moving  and  that 
T  is  moving  slowly  compared  to  the  measurements,  the  temporal  behavior  of 
the  received  signal  is  merely  that  induced  by  f  (t).  Specifically,  with  these 
assumptions  we  have 


b(0,  <j»,  t)  =  b(0 ,  <j>)  f(t  -  *) 


(4.  1) 


where  R  is  the  path  length  T5H.  This  then  represents  the  effect  at  R  at  time 
t  of  the  source  at  T  radiating  onto  dS.  The  term  b(0,  <j>)  represents  the  ran¬ 
dom  nature  of  the  signal  at  R  at  time  t. 

Consider  now  some  point  r  which  lies  in  the  receiving  aperture 
as  shown  in  Fig.  32.  Now  the  total  reflecting  surface  can  be  viewed  as  the 
source  of  many  waves  incident  on  £.  Specifically,  the  surface  dS  generates 
a  wave  incident  from  direction  k  as  defined  by  the  vector  from  it  to  dS.  Thus, 
we  can  consider  each  of  the  contributors  from  this  surface  to  point  r^  as  the 
contributor  of  a  plane  wave  with  an  amplitude  b(0,  4>)  or  equivalently  b(k) 
where  k  is  uniquely  defined  by  0  and  <$>.  This  allows  us  to  write  for  an  iso¬ 
tropic  (omnidirectional)  receiver  at  r  the  observed  signal  as 
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4} 


s(r,  t) 


■fh 


b(k)  exp(j  k  •  r)  f(t  -  — )  dk 
“  “  C  ™ 


(4.  2) 


For  simplicity,  let  us  assume  that  the  variation  of  f(t)  is  small  so  it  can  be 
neglected.  Thus  we  consider  s(r)  where 


b(k)  exp(j  k  *  r) 


(4.  3) 


Here  b(k)  equals  b(0,  <j>)  and  k  .  r  is  merely  the  projection  of  the  incident 
direction  vector  on  r. 

We  now  define  the  properties  b(k)  should  have  to  be  consistent  with 
properties  of  s(r).  Note  that  s(r)  is  a  random  field  in  r  generated  by  b(k). 
The  first  assumption  on  b(k)  is  that  it  is  a  zero  mean  Gaussian  random  field. 
Thus 


E[b(k)]=  0 


(4.4) 


This  assumption  is  consistent  with  the  roughness  assumption  leading  to  the 
definition  of  diffuse  multipath.  Since  it  is  Gaussian,  all  we  need  for  a  com¬ 
plete  statistical  description  is  the  covariance.  That  is,  we  need 


E[s(r)  sV)]  =  K(r,  r>)  . 


(4.5) 


This  is  written  as; 


Ks(r,  r1)  = 


r*. 


Ep>(k)  b  (k1)]  exp(j  k  •  r  -  j  k '  •  £')  dk  dk' 

(4.6) 
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It  is  convenient  at  this  point  to  assume  that  s(r)  is  a  homogeneous 
random  field  over  the  aperture  so  that; 

Ks(r,  r')  =  Ks(r  -  r')  .  (4.7) 

This  implies  (see  Yaglom)  that; 

E[b(k)  b*(k')]  =  K(k)  6(k  -  k1)  (4.  8) 

where  6(k  -  k')  is  a  two  dimensional  implse  function.  The  function  K(k)  or 
K(0,  <())  is  called  the  channel  spread  function  and  it  is  a  measure  of  the  energy 
coming  from  direction  (0,  <j>)  to  the  point  r. 

It  is  the  channel  spread  function  that  defines  the  angle  estimation  per¬ 
formance  of  antennas  (see  McGarty  [1])  and  thus  plays  a  dominant  role  in  deter¬ 
mining  the  effects  of  random  multipath  phenomena.  It  is  the  evaluation  of  this 
function  from  first  principles  and  the  use  of  it  in  system  performance  evalua¬ 
tion  which  will  interest  us  ir  as  section. 

4.  2  THE  CHANNEL  SPREAD  FUNCTION 

To  evaluate  the  effect  of  diffuse  radiation  on  azimuth  estimation  per¬ 
formance,  it  is  necessary  therefore  to  determine  the  channel  spread  function. 
To  obtain  this  function,  we  must  first  evaluate  the  time  behavior  of  the  irra¬ 
diated  surface  and  also  the  radiation  pattern  of  incremental  portions  the 
surface.  Once  these  two  relationships  are  obtained,  the  spread  function  can 
be  obtained  directly. 
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Consider  a  point  source  located  at  position  (0,  r,  z^)  on  the  flat 
earth  model  as  shown  in  Fig.  33.  A  receiver  is  located  at  (0,  0,  Zj).  The 
source  emits  a  spherical  wave  which  begins  at  time  tg  and  is  received  by  the 
receiver  at  time  tj  where 


t 


1 


(4.9) 


where  dg  is  the  distance  between  the  source  and  the  receiver  and  is  given  by; 


d 


0 


+  (z. 


(4.  10) 


The  power  density  of  the  wave  received  at  the  receiver  is  Pq/ (4  trdg)  where  P 
is  the  power  of  the  source. 

The  source  also  irradiates  a  portion  of  the  surface  of  the  earth  and 
this  surface  scatters  the  radiation  in  all  directions.  The  fraction  of  the  inci¬ 
dent  power  scattered  to  the  receiver  from  the  ith  incremental  area  is  given 
by  (Tj,  the  scattering  cross-section  per  unit  area.  Assume  that  there  are  N 
such  areas.  Then,  the  total  amount  of  power  per  square  meter  received  at 
the  receiver  due  to  the  N  discrete  incoherently  radiating  surface  areas  is, 

S^,  where; 


z 


Fig.  33.  Geometry  of  direct  path, 
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where  is  the  area  of  the  ith  scattering  section,  is  the  distance  from  the 
receiver  to  the  ith  scatterer,  and  is  the  distance  from  the  transmitter  to 
the  ith  area. 

Now,  at  any  one  instant  of  time,  not  all  of  the  ground  is  illuminated 
by  the  source.  For  example,  at  t  =  t^  the  wave  has  just  begun  to  propagate. 
The  minimum  amount  of  time  it  takes  for  a  wave  to  propagate  from  the  source 
to  the  receiver  via  a  reflection  from  the  earth's  surface  is, 


fr2+(Zl  +  z2)2 


(4.  12) 


where  c  is  the  speed  of  propagation.  At  this  instant,  the  area  of  the  scatter¬ 
ing  plane  illuminating  the  receiver  is  only  a  point.  For  t  >  t  ,  we  can  find 

s 

the  area  by  finding  the  equations  of  all  (x,  y)  such  that  the  distance  Rj  +  R^ 
equals  ct.  The  resulting  area  is  that  area  from  which  the  receiver  obtains 
power  at  time  t.  This  follows  directly  from  the  fact  that  Rj  is  the  distance 
from  the  scattering  point  to  the  receiver  and  R^  the  distance  from  the  scatter¬ 
ing  point  to  the  source. 

From  the  geometry  of  Fig.  34,  it  is  clear  that 


D2  2  ,  2,  2 

Rj  =  x  +  y  +  z  ^ 


(4.  13) 


R2  =  x2  +  (r  -  y)2  +  z2 


<4.  14) 


*•  f 
4) 


Then,  if  we  define  for  each  time,  t,  a  distance  p  as; 


p  =  ct 


(4.  15} 


we  can  show  that  the  area  irradiated  which  gives  rise  to  signals  received  at 
the  receiver  at  time  t  >  t  is  bounded  by  the  ellipse  given  by  (see  Appendix  I) 

S 


where 


2.  .2  ,  ,  2  2  2 

a  (y  -  y0)  +  b  x  =  rQ 


(4.  16) 


2  ,2  2. 

a  =  r(p  -  r  ) 


b2  =  4p2 


I  2  2,2  2\ 

.  I"  +Z1  "IP 


"2 — r 

2(p  -  O 


(4.  17) 

(4.  18) 

(4.  19) 


and 


2  2 
r0=  P 


2\2  A  2.  2  2. 

’  1/  '  4z2(P  '  r  ) 


T~2 - 2 7 

(p  -  r  ) 


(4.  20) 


Thus,  for  any  time  t,  the  surface  from  which  radiation  is  received  is  bounded 
by  the  ellipse  given  by  (4.  16).  We  shall  call  the  area  of  the  elliptical  surface 
A(t).  It  is  easily  shown  that  as  t  -*■  *>,  the  elipse  turns  into  a  circular  region. 
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Now,  if  we  assume  that  the  receiver  is  omnidirectional,  then  we  can 
calculate  the  total  power  density  at  the  receiver.  This  can  be  obtained  directly 
from  (4.  11)  as  a  limiting  case.  Namely,  if  <r(x,  y)  is  the  fraction  of  incident 
power  scattered  from  (x,  y)  into  the  receiver,  then  the  total  power  density  at 
the  receiver  at  time  t  is  given  by  S(t)  where 


•S(t)  - 


(4tr 


c r(x, 


A(t) 


y)  dx  dy 


(4.  21) 


where  A(t)  is  the  bounding  ellipsoid. 

There  is  an  important  interpretation  of  £q.  (4.  21)  worth  noting.  Con¬ 
sider  the  special  case  of  a  reflector  which  scatters  the  incident  radiation  uni¬ 
formly  in  all  directions.  For  this  case  we  let  cr(x,  y)  equal  <Tq  and  obtain  for 
S(t); 


S(t)  = 


dx  dv 


(4.  22) 


S(t),  now  in  Eq.  (4.  22),  is  the  power  per  square  meter  at  the  receiver  at  time 
t.  As  t  increases,  the  ellipse  grows;  however,  the  inverse  di  tance  squared 
values  weight  contributions  less  and  less.  It  can  be  shown  that  (4.  22)  approaches 
a  limiting  value  at  t  -*■  <*>.  An  exact  calculation  is  shown  in  Fig.  35.  Here, 
we  have  plotted  S(t)  versus  time  for  r  =  15  km,  z^  =  10  m  and  =  3  km. 

Note  how  rapidly  S (t)  approaches  almost  a  si.  ady- state  value.  What  this 
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Fig.  35.  Area,  S(t),  vs  time  (z,  =  10  m,  z,  =  3  km,  R  =  15  km) 


implies  is  that  there  will  be  little  time  fluctuation  in  the  diffuse  field  due  to 
this  effect. 

The  response  S(t)  due  to  pulses  is  shown  in  Figs.  36  and  37.  Here, 
the  source  emits  a  1  (isec  and  1.6  (J.sec  pulse  respectively.  The  correspond¬ 
ing  power  density  is  plotted.  Again,  note  the  sharp  rise  and  cutoff  values  for 
the  diffuse  radiation.  The  rise  time  is  about  0.  5  psec  while  the  decay  time  is 
about  0.  25  |j.sec. 

To  complete  the  analysis,  the  scattering  cross  section  p?r  unit  area 
<r(x,  y)  must  be  determined.  To  do  this,  we  shall  follow  Kodis  and  Barrick 
and  obtain  the  function  cr(0  ,  <}>  ;  $.).  This  is  the  fraction  of  the  incident  power 
incident  on  a  surface  S  from  and  angle  ^  and  scattered  into  direction  (6g,  4>g). 
These  directions  are  shown  in  Fig.  38.  Here,  it  is  assumed  that  the  incident 
field  is  a  plane  wave,  making  an  angle  <Jk  with  the  z  axis  and  having  a  wave 
vector  lying  in  the  y-z  plane. 

To  obtain  this  function,  we  write  the  scattered  field  as 


E  (r)  = 
— sc  — 


juji 


r(r,  r')  J(r')  dS 


(4.  23) 


where  T(r,  r1)  is  the  dyadic  Green's  function  given  by  (see  Silver,  p.  132) 

£(r,  r1)  =  (l  +  -j  V  vj  (exp[j  k|r  -  r'|  ]/(4ir|  r  -  r 1 1 ))  (4.  24) 


and  J(r’)  is  the  surface  current  density  and  j  is  J-l. 
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Fig.  36.  Differential  area  vs  time  (S(t  +  T)  -  S(t) ; 
T  =  1  Msec;  z,  =  10  m,  z,  -  3  km.  R  =  15  km). 
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Fig.  37.  Differential  area  vs  time  (S(t  +  T)  -  S(t); 
T  =  1.6  jusec;  z ^  =  10  m,  z^  =  3  km,  R  =  15  km). 


INCIDENT  WAVE 


Using  the  Kirchoff  approximation,  Kodis  shows  that  the  field  scattered 


in  direction  k  due  to  a  source  from  k.  is  given  by 


E  (k.,  k  )=J^  ejkrf  dS'(I  - 
— sc— 1  —s'  2irr  I  - 


kg  kg)  (h  X  ft) 


exp [j(k^  -  1^)  •  r'j 


(4.  25) 


where 


h  =  k.  x  e 


(4.  26) 


and  his  the  unit  normal  directed  into  the  scattering  plane  at  every  point  and 


eis  the  polarization  vector  of  the  incident  field.  Using  the  method  of  station¬ 


ary  phase  and  then  averaging,  Kodis  shows  that  the  scattering  cross  section 


is  given  by 


o-(es.  <|>8;  4>.)  =  irnA<  |rj  r2|>|R(?)|‘ 


(4.  27) 


where  nA  is  the  average  number  of  specular  points  per  unit  area,<  jr^  r^j  > 
the  average  absolute  value  of  product  of  principal  radii  of  curvature,  and  R(?) 


the  reflection  coefficient.  The  angle  ?  is  the  local  angle  of  incidence  at  the 


specular  point  and  is  given  by  the  relationship; 


cos?  =  ^  Jl  -  sin  sin  cj>8  cos  0^  +  cos  4>^  cos  4>s  4  (4.  28) 


y  < 


1 


% 


s 


Barrick  has  shown  that 


nA  = 


7.  255 

2  '2 
tt  2<r. 


(4.  29) 


with 


2 

s 


=  4 


2 

*h 


(4.  30) 


where  <r^  is  the  correlation  length  of  the  surface  and  is  the  variance  of  the 
surface  height,  tan  y  is  given  by 


tan  y  = 


I.  2  2 

J sin  cj>^  -  2  sin<j)  sin<j>g  cos0g  +  sin  4^ 
cos^TTToi^T 


(4.  3  ’ ) 


The  curvature  has  also  been  shown  by  Barrick  to  be 


<r,  4 

<  |rj  r2(>  =  0*  1387  it  — ^  sec  y  , 


(4.32) 


Finally,  for  vertical  polarization  over  a  perfectly  conducting  plane,  we  have 
for  the  reflection  coefficient 


-  sin$,  sin^  sin4  6  +  a_  a 

R(C)  = - - - U - £ - — 

4  sin4?  cos4? 


(4.33) 
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where 


a  -  =  cos<j>.  sinA  +  sin<j>.  cos<j>  cos0 

u  X  S  X  S  £ 


(4.  34) 


a.  =  sin<|>.  cosifi  +  cos4>.  sin<}>  cosG 
j  IS  X  s  s 


(4.  35) 


If  we  define  a  constant  C  as  the  product  of  0.  1378  and  7.  255  then  we 
can  combine  (4.  28)  -  (4.  35)  to  yield; 


o-fV  4>s;  <t>i)  =  C  —  y  exp 


tan^y 

'  2 
s  . 


R(5) 


(4.  36) 


This  is  the  resulting  scattering  cross-section.  We  can  observe  its  behavior 

easily  in  this  form.  By  using  (4.  31)  we  first  see  that  for  positive  <J>.  and  <j>  , 

X  s 

tan  v  increases  as  the  azimuth  0  increases  or  decreases  from  zero.  This 

implies  that  the  angle  y  is  increasing  as  an  absolute  value  of  the  azimuth. 

Thus  for  fixed  <f>^  and  the  azimuth  behavior  of  the  scattering  cross-section 

is  dominated  by  y  which  in  turn  is  reflected  in  the  two  terms;  exp(-tan  y/a  ) 

2 

and  sec  y.  The  former  term  decreases  with  increasing  y,  the  rate  of  increase 

2  4 

depending  on  the  value  of  s  .  The  latter  term,  increases  with  sec  y,  inde- 

2 

pendent  of  s  .  Now  two  distinct  regions  are  possible.  The  first  is  for  large 

4 

s;  In  that  case,  the  sec  y  behavior  dominates  a  scattering  cross-section 
which  increases  with  azimuth.  Large  s  implies  thatcr^,  the  correlation 
length,  is  small  compared  to  the  surface  roughness.  This  can  be  called  a 
very  rough  surface.  For  example,  if  is  10  m,  then  the  value  of  cr^  needed 


to  insure  this  condition  is  on  the  order  of  10  m.  A  second  effect  to  note  is 
that  the  amplitude  about  zero  azimuth  depends  inversely  on  s^.  Thus,  a  very 

rcugh  surface  scatters  less  in  the  forward  direction.  The  second  region  is 

2  2  2 
that  in  which  s  is  small,  i.  e.  ,  cr^  »  cr^.  so  that  exp(-tan  y/s  ;  dominates 

the  profile  of  the  scattering  cross-section.  Figs.  39-42  depict  results  for 
varying  cases.  Note  that  in  Fig.  40  we  have  almost  a  specular  scatter  sur¬ 
face  with  minimum  spread.  This  results  from  a  very  large  value  of  the 
surface  correlation  length. 

The  assumptions  that  must  be  made  concerning  the  surface  in  order 
that  the  analysis  leading  to  (4.  34)  be  consistent  are  (see  Barrick): 

1.  The  radius  of  curvature,  p,  everywhere  on  the  reflecting 
surface  be  much  greater  than  the  wavelength.  Thus  p  >>  \q. 

2.  Multiple  scattering  effects  can  be  neglected. 

j.  The  mean  square  surface  height  is  much  greater  than  a 
wavelength,  that  is  cos  ^  >> 

Condition  (1)  and  (3)  are  direct  limitations  on  the  types  of  surfaces. 
Specifically,  (1)  requires  that  the  surfaces  be  very  rounded  since  X.q  at  1090 
Mhz  is  about  one  foot.  Furthermore,  (3)  states  that  unless  h  is  very  large, 
the  approximations  breakdown  at  grazing  angles;  <j>j  -*■  tr/2.  That  is  most  sur¬ 
faces  "look"  like  specular  reflectors  at  grazing  angles  such  that  4>^  a>  it/2* 

Further  limitations  concerning  the  surface  correlation  coefficient 
are  discussed  in  Barrick  [2].  He  notes  that  the  correlation  must  be  quadratic 
about  the  origin.  Specifically,  the  result  in  (4.  36)  assumes  that  the  surface 
height  z  is  a  Gaussian  random  variable  with  a  correlation,  E[z(x,  y)  z(x',  y ' )], 
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Fig.  40,  Scattering  cross-section  vs  azimuth  (0.  =  <f>  -  88°,  s  =  1/2  ). 
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having  the  form,  cxp{-  a  p  ),  where  p  is  (x  -  x'y  +  (y  -  y’)  .  This  form 
ensures  power  conservation  and  correct  behavior  as  the  ray  optics  limit  is 
approached. 

^  We  now  want  to  evaluate  the  spread  function  K(9,  <|>)  which  represents 

the  power  per  unit  area  per  square  radian  incident  at  the  receiver  coming 
from  directions  $  in  elevation  and  6  in  azimuth.  To  do  this,  consider  an  incre¬ 
mental  area  AS'  as  shown  in  Fig.  43.  Draw  a  line  from  the  origin  at  an  angle 
9  in  the  x,  y  plane,  bounding  dS'  on  one  side  and  a  line  at  angle  9  +  A 9  bound¬ 
ing  AS'  on  the  other  side.  Now  <j>  is  given  by 


4>  =  tan 


-1 


j~Z  ~  2 

V*  l  y 


(4.37) 


The  loci  of  constant  for  constant  are  circles  in  the  x,  y  plane. 
Circles  for  $  and  4>  +  A<p  are  shown  also  in  Fig.  43.  The  cross  hatched  area 
in  this  figure  represents  the  region  of  the  scattering  plane  from  which  the  source 
receives  radiation  coming  from  the  solid  angle  $  +  A<j>,  and  0,  0  +  A  9. 

The  power  per  unit  area  received  from  this  area  is  K(8,  (f>)  A 9  A<j>.  Also, 

AS  is  directly  defined  by  the  bounding  region. 

For  any  point  in  the  area  from  Fig.  44,  we  see  that 


(4.  38) 


The  scattered  angle  4>  is  given  by  (4.  37).  The  angle  8,  the  angle  used  for 

8 

the  scattering  cross-section  is  given  by 


~  "r 


e*  =  e  +  4> 


(4.  39) 


where 


<p  -  tan 


-If  x 
r  -  y 


(4.  40) 


Thus,  we  have  for  any  0,  <j>  pair 


K(0,  <j>)  =  lim  (4ir)  j 
A6  -  0  Q 

A4>  -  0  & 


0  C  I  cr(0  +  <p,  <j>;  4>.)  dx  dy 

?JJ  44 


Ad  Ac}> 


(4.41) 


where  $  is  the  intersection  of  the  region  formed  as  in  Fig.  43  and  the  area  of 
illumination  at  time  t.  Note  that  x,  y  and  4>^  are  functions  of  0  and  <}>  so  that 
by  using  the  Jacobian  of  the  transformation  some  simplification  of  (4.  41)  may 
result.  However,  another  simplifying  tool  can  be  used  if  we  introduce  variable 


Pq  where 


I  2  7  2 

p0=  yjx  +  y 


But,  from  (4.  37)  and  Fig.  44  we  note  that 


(4.  42) 


pQ  *  tan  cj> 


(4.  43) 


i 


Here  pg  is  the  radius  of  a  circular  coordinate  system  in  the  x,  y  plane.  Then 


x  =  pg  sin8  =  z  j  tan4>  sin0  (4.  44) 


y  =  Pg  cos0  =  Zj  tan4>  cosfi  .  (4.45) 


With  these  substitutions  we  can  write  ip  as  in  (4.  40); 


ip  =  tan 


-1 


Zj  tan$  sin9 


r  -  Zj  tan<j>  cos  6 


(4.46) 


and 

Rj  =  Zj(l  +  tan24>)  (4.47) 

R2  =  z2  tan2(jj  +  r2  -  2r  z^  tan<j>  cos0  +  z^  .  (4.  40) 

Then  using  the  transformation  from  x,  y  into  pg,  8  the  integral  in  (4.  41)  can 
be  written  using 


dx  dy  =  p0  dpQ  d0 


(4.  49) 


as; 
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(4.  50) 


K(6,  <j>)  =  lira  — - 

A0  -  0  (4ir) 
A<)>  “*•  0 


0  II  <r(B  +  ip(d,  <*>),  <J>;  <^(0, 

r)2  n 


4»))p0  dP0  de 


A  0  A4> 


where  we  have  assumed  that  A(t)  is  large  enough  (t  >>  t^)  and  have  explicity 
stated,^,  Rj  and  R2  in  terms  of  0  and  tp.  Now  using  (4.  49)  and  (4.  43)  we  can 
take  the  limit  to  yield; 

Prt  cr(0  +  4>),  4>;  4>.(0,  <j>))  zf  sec 2  <j>  tan4> 

K{0,  <>)  =  - - - - i - - - - - (4.51) 

m  Rj(6»  <(>)  R^Q.  4>) 

as  the  resulting  channel  spread  function  for  this  diffuse  multipath  field. 

An  evaluation  of  the  channel  spread  function  has  been  performed  and 
the  results  are  shown  in  Figs.  45  to  48.  This  is  for  a  time  large  enough  so 
that  all  the  surface  is  irradiated.  The  ground  antenna  is  assumed  to  be  at 


20  m  and  the  aircraft  at  1  km.  The  value  of  s  varies  as  does  the  range  of  the 


2 

aircraft.  In  these  figures,  we  have  normalized  (4.  51)  by.P0/4ird0  and  have 
plotted  10  logjQ  of  that  ratio. 

In  Fig.  45,  we  have  a  fairly  smooth  surface  (s  =  0.  5).  We  note  the 
bright  spot  near  the  horizon  and  a  rapid  decrease  as  6  goes  from  broadside. 

A  similar  plot  is  shown  in  F ' hut  now  the  aircraft  has  quintupled  its 
range.  We  note  the  narrowing  of  the  peaks.  This  narrow  area  is  the  glisten¬ 
ing  surface  mentioned  in  Beckmann  and  Spichizzino.  The  remaining  two  plots 
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show  how  sidelobes  appear  as  the  surface  becomes  more  rough.  This  is  a 
result  of  the  increase  in  the  radius  of  curvature  term  as  well  as  a  slow 
decrease  in  the  average  number  cf  scatters  as  appear  in  (4.  51). 

With  K(6,  <j>),  the  azimuth  spread  for  a  given  antenna  vertical  gain 
can  be  evaluated.  Let  G(<j>)  be  the  vertical  gain  of  the  antenna;  then,  azimuth 
spread  K(0)  is  given  by 


K(0)  = 


ir 

'1 

K{0,  <j>)  44#  . 


(4.52) 


As  before,  this  is  a  function  of  cr^,  cr^,  z,,  z^,  r,  and  time.  A  sample  plot 
for  r  =  30  km  using  the  vertical  array  proposed  for  the  DABS  antenna  is 
shown  in  Fig.  49.  Note  the  sharp  central  peak  and  the  drop  of  6  orders  of 
magnitude  as  the  azimuth  increases.  This  profile  is  called  the  glistening 
surface  by  Beckmann  and  Spizzichino. 

It  is  also  possible  to  use  the  gain  in  azimuth.  Thus,  if  G(0,  c)>)  is  tie 
gain  of  an  antenna  in  elevation  and  azimuth  then,  the  quantity 


IT 

1 


ir 


G{0,  <)))  K(0,  4>)  d0  d<js 


n  = 


ir 

"2 


/-IT 


4ir  d. 


(4.  53) 
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Fig.  49.  Plot  of  channel  spread  function  vs  azimuth  for  two  correlation  lengths 


The  ambiguity  function  of  array  is  defined  as^$  (p)  and  is  given  by 


B  fa')  =  m(p  +  H-') p 


where  m(p)  is  the  N  x  1  complex  vector 


m(p)  - 


exp(-j2«d^A0)‘ 


_exT>(-j2irdNnA0)J 


(4.  ST) 


(4.  58) 


Here,  d^  is  the  displacement  from  the  origin  of  the  ith  dipole  receiver  (for  a 
dipole  array),  is  the  wavelength,  and  N  is  the  number  of  dipoles. 

The  total  signal  energy  is  given  by 

<rf  =  cr2  E  (4.  59) 

8  D  8 

2 

where  Efl  is  the  energy  of  the  transmitted  pulse  and  represents  channel 
attenuation  effects.  The  total  diffuse  multipath  signal  level  is 

c r2  E  <r.2  (4.  60) 

1  si 
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since  we  assume  the  signal  due  to  diffuse  multipath  has  been  generated  by 
the  source.  The  noise  is  assumed  to  be  a  zero  mean  Gaussian  white  noise 
process,  independent  from  element  to  element  in  the  array  and  has  a  variance 
given  by 

(4.61) 

it  becomes  convenient  to  define  two  ratios,  the  multi  - 
and  the  noise -to -signal  ratio,  y  These  are 

(4.  62) 

(4.  63) 

Now  it  has  been  shown  elsewere  that  the  Cramer-Rao  bound  for  estimating  ji 
is  given  by  (McGarty  [1  ]); 

var(|x  -  pi)  >  -  j  (4.  64) 

where 


2  _  N0 
Si  *  T 


For  analytical  purposes 
path -to- signal  ratio,  ?, 


2 

<r 

Y  =  n 
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(4.70) 


0*0  =  -tan2  6  +  tan2  0 


1  +  2 


cot2  9 
sin2  0 


2 

<r 


where  0  is  the  pointing  angle.  Note  that  0  =  90°  implies  broadside. 

Now  the  evaluation  of  (4.  65)  can  be  facilitated  by  using  Fourier  trans¬ 
form  relatione liips.  This  has  been  discussed  elsewhere  but  the  evaluation  of 
the  performance  curves  in  Fig.  51  have  been  evaluated  using  this  technique. 
(See  McGarty  [1], ) 

We  can  express  the  performance  as  a  function  of  the  beamwidth  of 
the  array  as  follows.  Recall  that  for  a  linear  array,  we  have  for  the  ambiguity 
pattern 


8w  = 


sin 


Nd 

~T 


"7"<r . 

Sin  ry  ji 


(4.71) 


where  d  equals  Zirdp/X-Q,  d^  being  the  interdipole  spacing  and  \q  the  wavelength. 
The  first  null  in  this  pattern  occurs  at  where 

•  (4.72) 

The  center  of  the  beam  is  at  p  =  0  so  that  the  beamwidth  between  nulls 
is  given  by 


(4.73) 
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From  McGarty  [1],  we  have  the  performance  bound  in  the  absence 
of  multipath  which  yields  the  standard  deviation  of  the  azimuth  estimate,  o-0 
in  terms  of  the  signal -to -noise  ratio  and  the  number  of  elements.  It  is 


crfi(  radians)  =  — - — , - 

\'vr  (hr  -  N) 


(4.74) 


where  or  is  the  signal -to-nois©  ratio.  Now,  using  (4.73)  for  N  we  have 


(4.75) 


0  BJ 


k0  sluoB, 


Assuming  0^  <<  1,  we  can  approximate  this  by 


2/0 

1 


where  <Tq  is  in  radians  as  is  0^. 


(4.76) 


That  is,  in  the  absence  of  multipath  the  standard  deviation  of  the 
angle  estimate  is  proportional  to  the  three -halves  power  of  the  beamwidth. 
This  curve,  along  with  the  curve  including  multipath,  is  shown  in  Fig.  51. 

Also,  if  the  multipath  effect  is  small  (i.  e. ,  ?  >>  1),  then  the  error 
with  multipath  can  be  written  as; 


$  V# 
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=  o-g,  no  multipath 


1  +  S 


K(n')§£&  (ji  -  jx«)  d(J.*] 


(4.77) 


That  is,  the  variance  is  increased  by  a  factor  which  depends  on  the  multi  - 
path  to  signal  ratio,  5,  and  the  nature  of  the  spread  channel,  K(|jl). 

Plotted  in  Fig.  51  is  <Tg  versus  0g,  the  beamwidth  in  degrees  for  dQ 
equal  0.  152  m,  a  frequency  of  1090  MHz,  a  signal -to -noise  ratio  of  25  dB, 
and  the  interference -to -signal  ratio  of  -20  dB.  We  use  the  spread  function 
for  Zj  equal  to  10  m,  z 2  equal  to  2  km.  range  equal  to  20  km.  cr^  was  10  m, 
and  <r^  ,  200  m  (as  shown  in  Fig.  49).  Note  the  behavior  of  the  variance.  As 
0g  decreases,  c Tq  decreases  almost  linearly.  The  bound  actually  is  dominated 
by  S/N  behavior  and  is  only  slightly  affected  by  the  diffuse  multipath.  This 
implies  that  for  an  array,  with  an  8\q  vertical  aperture  and  a  26\.^  horizontal 
aperture  diffuse  multipath  will  not  be  a  significant  problem. 

4. 4  CONCLUSIONS 

Diffuse  multipath  can  be  modelled  as  a  random  phenomenon  which 
results  from  the  rough  and  random  nature  of  the  scattering  surface.  In  this 
section  we  consider  the  effects  of  such  randomness  in  detail.  We  first  reviewed 
the  rather  large  body  of  work  that  already  exists  on  scattering  from  random 
surfaces  and  discussed  the  shortcomings  or  usefulness  of  several  of  the  models 
that  have  been  proposed. 
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A  specific  model  for  diffuse  scattering  was  then  proposed.  First 
an  analysis  of  the  area  irradiated  as  a  function  of  time  was  performed.  The 
result  of  this  analysis  was  that  for  psec  type  signals  the  effects  of  finite  time 
propagation  could  be  neglected.  Then  using  the  Kodis-Barrick  model  the  scatter¬ 
ing  cross-section  per  unit  area  was  obtained  for  a  Gaussian  surface  with  a 
Gaussian  correlation  function.  The  resulting  scattering  nroes -section  cr C0 s  , 

4>g5  <J>^)  was  evaluated  in  terms  of  the  mean  square  surface  height,  cr^,  the 
correlations  length  of  the  surface,  o^,  and  the  geometrical  parameters  of  the 
point  of  scattering. 

It  has  been  shown  that  to  evaluate  the  effects  of  random  signals  coming 
from  a  continuum  of  directions,  as  is  the  case  for  diffuse  multipath,  the  channel 
spread  function,  K(8,  <j>),  must  be  obtained.  This  function  represents  the  frac¬ 
tion  of  diffuse  energy  received  from  a  given  azimuth,  0,  and  elevation,  cj>. 

Using  the  Kodis-Barrick  scattering  cross-section  and  the  geometry  of  the  ground- 
to-air  surveillance  link  the  spread  function  was  obtained. 

This  section  concluded  with  an  evaluation  of  the  effects  of  diffuse  multi - 
path  on  azimuth  eatimation.  Several  different  types  of  surfaces  were  considered 
and  a  simplified  closed  form  analytical  result  for  azimuth  errors  was  obtained. 

As  a  result  of  this  study  the  following  specific  conclusions  were 

reached. 

1.  The  Kodis-Barrick  model  for  diffuse  scattering  cross-sections 
is  most  appropriate  for  performance  evaluations. 

Time  dispersion  effects  are  not  important  for  MHz  data  rate 
transmissions  with  respect  to  diffuse  multipath  effects. 


2. 


The  channel  spread  function  for  cr^  <  <r ^  shows  a  distinct 
glistening  surface  effect  that  has  been  noted  by  other  investi¬ 
gators. 

The  errors  on  azimuth  estimation  due  to  diffuse  multipath  are 
negligible  and  orders  of  magnitude  smaller  than  those  observed 
for  specular  multipath.  Thus  diffuse  multipath  presents  no 
serious  limitations  in  direction  finding  capabilities  in  DABS. 
The  effects  of  earth  curvature  and  shadowing  can  be  included 
but  their  net  effect  is  of  second  order  and  discussed  in  the  next 
section. 


SECTION  FIVE 
OTHER  EFFECTS 


In  the  past  two  sections,  we  have  discussed  the  effects  of  specular 
and  diffuse  multipath  and  have  presented  models  with  which  the  effects  of 
these  sources  of  interference  on  system  performance  could  be  determined. 
In  so  doing,  we  neglected  certain  factors  which  in  some  circumstances  can 
become  important.  It  is  the  purpose  of  this  section  to  discuss  these  factors 
and  place  them  in  their  proper  perspective  relative  to  the  general  model. 


5.  1  REFRACTION 

The  model  developed  in  this  report  for  specular  multipath  considered 
the  scattering  plane  to  be  flat.  However,  the  earth  is  not  Uat,  it  has  curvature. 
If  we  consider  any  scattering  point  x,  y  and  let  r ^  he  the  distance  from  the 
transmitter  to  x,  y  and  r ^  the  distance  from  the  receiver  to  that  point,  then 
Beckmann  and  Spizzichino  (pp.  222-224)  show  that  we  must  reduce  the  reflec¬ 
tion  coefficient  derived  in  Section  3  by  a  factor  D  where 


D  =  1  + 


2rl  r2 


Re(rl  +  r 2'  3in  Y 


1/2 


(5.1) 


where  Rg  is  the  effective  radius  of  the  earth  and  y  is  (ir/2)  -  4^  where  4^  is  as 

in  Section  4.  Then,  for  r, ,  r,  «  R  ,  this  is  close  to  unity  and  it  can  be 

Z  2 

neglected.  In  our  analysis,  we  had  also  included  the  attenuation  of  1/r^r^ 

which,  combined  with  D,  would  imply  for  large  R  that 

© 


( 


s 


i 
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(5.  2) 


rl  V Z  rl  r2 


For  the  angles  of  interest  (y  *  3  )  and  r^,  r^  <  10  km,  D  is  sufficiently  close 
to  one  to  be  neglected.  However,  for  a  more  general  analysis,  its  effect  must 
be  included. 

Tn  the  same  context  of  this  problem  is  the  evaluation  of  the  line- of  - 
sight.  Namely,  for  a  given  antenna  height  and  aircraft  height,  how  far  from 
the  antenna  must  the  aircraft  fly  to  fall  below  the  horizon.  If  z^  «  z^  then, 

R^,  the  range  to  the  horizon  can  be  shown  to  be 


Rh  “  JZz  2  RE 


(5.3) 


where  z^  is  the  altitude  of  the  aircraft  and  Rg  the  effective  radius  of  the  earth. 

Rjj,  equals  4/3  Rg  where  Rg  is  the  geometrical  radius  of  the  earth.  A  plot 
o  o 

of  R^  vs  z^  appears  in  Fig.  52.  As  seen  for  an  altitude  of  1  km,  the  range  to 
the  horizon  is  130  km. 

5.  2  SHADOWING 


A  second  effect  to  be  considered  is  that  of  shadowing.  Here,  we  are 
interested  in  the  decrease  in  area  irradiated  due  to  the  surface  irregularities 
themselves.  This  was  first  discussed  by  Beckmann  [3]  who  develop"'*  an 
approximate  expression  for  the  fractional  effective  area.  Due  to  an  error  in 
his  analysis  (see  Shaw),  Beckmannts  theory  underestimates  the  effects  of 
shadowing.  A  geometrical  approach  was  suggested  by  Brockelman  and  Hagfors 
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Fig.  52.  Plot  of  aircraft  altitude  vs  maximum  range  to  horizon. 


and  |(t)  is  the  surface  height  (e.  g.  ,  «(x,  y)).  Depending  on  the  surface  correla¬ 
tion,  (5,6)  and  (5.7)  are  directly  evaluated.  Wagner  shows,  that  for  the  range 
of  values  discussed  in  Section  4,  S(4>j,  c^)  is  about  0.  1  This  indicates  an  order 
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of  magnitude  decrease  in  the  values  of  the  field  intensity  evaluated  there. 

Since  there  appears  to  be  little  data  to  subtantiate  these  results,  we  did  not 
include  them  in  the  analysis. 

5.  3  SEA  SURFACE  SCATTERING 

The  scattering  of  electromagnetic  waves  from  a  sea  surface  is  much 
more  complex  than  that  from  land.  A  great  deal  of  both  the  theoretical 
and  experimental  work  has  been  done  in  this  area  but  none  as  of  yet  directly 
applicable  to  the  DABS  problem.  The  theoretical  work  by  Beard  [l]-[2], 

Beard,  et  al.  ,  Barrick  [1]- [3]  and  the  experimental  work  of  Bass,  et  al.  ,  [1]~ 
[2]  provide  some  insight  into  the  problem.  In  particular,  it  has  been  noted 
that  specular  type  reflections  can  arise  in  other  than  the  expected  specular 
direction.  Furthermore,  the  difference  between  diffuse  and  specular  reflec¬ 
tions  is  not  as  clear  as  that  developed  in  Sections  3  and  4. 

The  work  by  DeLorenzo  and  Cassedy  has  shown  how  one  might  analyze 
and  interpret  the  results.  Their  model  assumes  that  the  surface  of  the  ocean 
is  a  narrowband  zero  mean  Gaussian  random  process.  They  assume  that  the 
height  of  the  ocean,  z(x,  y),  depends  only  on  x  and  that  (see  Fig.  53) 


E[z(x)]  =  0 


(5.  10) 


Efz(Xj)  z(x2)]=  Kz(xj  -  x2)  . 

The  field,  z(x),  is  homogeneous. 


(5.11) 
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Using  a  Floquet  analysis  based  upon  Cassedy,  they  proceed  to  identify 
the  specular  nature  of  the  average  intensity  of  the  scattered  field.  They  do 
this  by  assuming  that  a  length  2L  of  the  ocean  is  illuminated  and  write  the 
normalized  scattered  field  as; 


P  = 


exp(jv  x  -  ji/  z(z))  dx 

X  s-t 


(5.12) 


where 


v  -  k(sin«j>_  -  8in4>  )  (5.  13) 

X  us 

s 

V 

a' 

v  =  k(co;v<|>Q  +  cos<j>g)  (5.  14) 


where  k  equals  2tt/X.q.  Clearly,  =  0  defines  the  classical  specular  direc¬ 
tion  since  <J>n  is  the  angle  of  incidence  and  $  the  scattering  angle, 
u  s 

By  assuming  that  Kz(r)  equals 


kz(t)  =  exp(-  t2/tq)  cos(2Trfgr) 


(5.  15) 


2 

where  is  the  variance  of  the  surface  height,  Tq  the  correlation  length  and 
f  the  spatial  periodicity  of  the  waves,  they  find  that  for  the  case  Tq  »  L; 

2  2  V"'  2  2  2  8ln2^x +  n2irfs)L 

E[pp^]  =  exp(-  v\  <r2 )  .  2^  ln<vzv\J  4L  -J— 


ns: -00 


[(vx+  n2irfg)L] 


(5.  16) 
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This  expression  shows  the  scattered  power  to  be  strongly  peaked  at 

v  +  2irnf  =  0  (5.  17) 

Obviously,  v  ~  0  is  one  of  these  values  but  so  too  are  other  directions. 
DeLorenzo  and  and  Cassedy  show  data  which  substantiates  their  model.  Other 
variations  of  the  surface  parameters  introduce  diffuse-like  scattering  along 
with  multiple  specular  scattering  directions. 

Thus,  the  analysis  of  the  effects  of  multipath  over  the  sea  appears  to 
be  quite  complex  and  still  an  open  area  on  many  issues. 

5.  4  COMPOSITE  SCATTERING 

The  fourth  and  final  extension  would  be  the  scattering  from  composite 
surfaces  as  discussed  by  Beckmann  [1].  Basically,  a  composite  rough  surface 
is  one  having  several  cr^  and  cr ^  Beckmann  shows  that  the  total  energy  from 
such  a  surface  is  the  sum  from  all  of  its  composites.  To  do  this  in  practice 
requires  a  better  understanding  of  the  surface  topography.  The  studies  by 
Fung  and  Moore  and  Hayre  and  Moore  do  not  seem  to  indicate  however,  that 
such  composite  behavior  is  common. 

5.  5  CONCLUSIONS 

Those  other  effects  that  were  not  considered  in  the  analysis  of  specular 
and  diffuse  multipath  in  Sections  3  and  4  respectively  represent  only  second 
order  perturbations  on  the  results.  The  earths  curvature  giving  rise  to  diver¬ 
gence  alters  the  result  only  slightly.  The  only  serious  limitation  of  having  a 
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curved  earth  is  that  of  line  of  sight  limitations.  These  can  be  easily  included 
as  a  cutoff  at  maximum  range. 

The  effects  of  shadowing  are  also  only  second  order  and  only  tend  to 
further  reduce  the  effects  of  diffuse  multipath.  However,  sea  surface  scatter¬ 
ing  does  present  a  phenomenon  for  which  little  data  in  the  L  band  regions  exist. 
This  effect  will  undoubtedly  be  important  at  airports  which  are  adjacent  to  the 
ocean.  To  a  first  order  the  sea  surface  can  be  considered  as  a  slightly  rough 
flat  plane  with  no  tilts.  This  would  indicate  that  azimuth  errors  would  be  small. 

There  are  two  other  effects  which  may  be  dominant  in  certain  locations 
and  for  which  at  present  little  is  known.  They  are; 

1.  Vegetation  Effects:  The  effect  of  trees  and  other  forms  of 
vegetation  must  be  considered  when  evaluating  the  reflection 
coefficient.  However,  trees  at  L  band  represent  a  random, 
finite  thickness,  lossy,  dielectric  scatterer.  At  present  no 
adequate  theory  is  available  to  analyze  this  effect.  DABSEF 
data  indicate  that  a  reduction  in  reflection  coefficient 
results. 

2.  Low  Angle  Effects;  When  the  aircraft  is  at  very  low  eleva¬ 
tion  angles,  theory  for  flat  surfaces  indicate  that  p  approaches 
unity  and  very  deep  nulls  occur.  However,  actual  topography 
is  not  perfectly  flat.  There  are  many  different,  albeit  large 
in  area,  surfaces  of  various  tilts.  Thus  as  the  elevation 
angle  gets  smaller  the  Fresnel  zone  gets  larger  and  encom¬ 
passes  many  independently  oriented  surfaces.  The  coherent 
sum  of  the  plane  waves  from  these  surfaces  tend  to  cancel  and 


the  resultant  deep  nulls  do  not  appear.  Again  no  adequate 
theory  exists  to  model  this  phenomenon  nor  are  adequate 
data  available  on  large  scale  correlation  effects. 


SECTION  SIX 


APPLICATIONS 

The  models  that  have  been  developed  for  both  diffuse  and  specular 
multipath  are  useful  in  evaluating  many  of  the  important  aspects  of  DABS.  In 
this  section  we  shall  discuss  five  areas  in  which  these  analytical  methods  find 
application  and  discuss  how  the  results  may  be  applied.  These  areas  have  in 
common  an  interest  in  how  the  signal  behaves  as  a  function  of  sensor  type, 
sensor  location  and  aircraft  location.  In  general,  the  resulting  analyses  and 
conclusions  will  be  based  solely  upon  the  specular  multipath  model,  for,  as 
we  have  shown  in  Section  4,  diffuse  multipath  is  a  second  order  effect. 

6.  1  PERFORMANCE 

There  are  two  aspects  of  DABS  performance  that  are  impacted  upon 
by  multipath.  The  first  is  the  communication  link  performance.  The  link 
performance  depends  upon  SNR  which  in  turn  can  be  seriously  degraded  when¬ 
ever  the  signal  fades  due  to  coherent  multipath  cancellation.  The  depths  of 
these  fades  can  be  evaluated  for  various  terrain  conditions  and  the  link  perfor¬ 
mance  obtained.  In  Section  3,  we  found  that  fading  depended  upon  geometry, 
surface  roughness  antenna  gains  and  sensor  and  aircraft  locations.  In  con¬ 
trast  to  a  simple  flat  earth,  single  plane  model,  the  model  developed  in  Sec¬ 
tion  3  was  based  upon  a  more  robust  set  of  conditions.  It  allows  fade  calcu¬ 
lations  and  performance  predictions  to  be  made  from  multiple  reflecting 
planes  of  finite  sise  and  varying  surface  interface  properties. 
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The  second  performance  issue  is  that  of  azimuth  estimation  perfor¬ 
mance.  As  was  shown  in  Section  3,  the  errors  in  monopulse  azimuth  estima¬ 
tion  are  dependent  upon  two  effects.  The  first  was  a  pure  SNR  fade  and  the 
second  was  a  coherent  cancelling  by  the  off-azimuth  multipath  plane.  Simi¬ 
larly,  ATCRBS  angle  estimating  schemes  have  also  been  analyzed  and  similar 
error  behavior  has  been  observed. 

6.  2  LOCALIZATION 

Section  2  was  basically  an  analytical  approach  to  specular  multipath 
in  that,  for  a  given  geometric?1  configuration,  the  resulting  output  signal  could 
be  formed.  The  inverse  of  this  technique,  going  from  an  observed  output  to 
a  given  geometry,  is  the  synthesis  portion  of  the  analysis  and  has  been  termed 
localization.  That  is,  by  observing  the  output  of  an  antenna  and  knowing  the 
model  it  is  possible  to  obtain  the  angles  a. .,  &  as  well  as  the  offset  location  of 
the  planes  giving  rise  to  that  signal.  In  such  a  fashion  we  can  localize  the 
source  of  each  specular  multipath  signal  and  identify  it  in  the  terrain. 

From  Section  3  we  know  that  the  specular  multipath  signal  power, 
sampled  in  time  and  normalized  by  the  intensity  of  the  direct  path,  is  given  by 

n  n 

s2  =  1  +  Yi  p 2  +  2^P  Pi  cos(6 
i=  1  i=  1 

n  n 

+  *YY*  Pi  cos (6^  -  . 

i=i  j=l 


(6.1) 


*r  »r  •  ■ 


Here  the  represents  path  effects,  range  effects  and  antenna  effects.  It  is 
not  the  reflection  coefficient  of  Section  3  but  a  conglomerate  of  all  effects  nor¬ 
malized  by  the  conglomerate  of  all  direct  path  effects.  The  factor  n  represents 
all  paths  except  the  direct  free  space  path.  The  phase  6^k  representes  the 
phase  difference  due  to  geometry  between  the  direct  path  and  the  kth  multipath. 
Note,  that  for  convenience,  we  have  neglected  the  Fresnel  reflection  coeffi¬ 
cient  phase  factor. 

Let  the  aircraft  be  at  a  position  (r  siny,  r  cosy,  z^)  and  the  receiver 
at  (0,  0,  Zj)  where  r  is  the  horizontal  ground  range.  Then  for  reflecting 
planes  with  small  range  and  cross  range  slopes  (less  than  about  +  10°)  it  can 
be  shown  that  6</'k  is 

2irf. 

S'/'k  =  -7-  +  \  (6-  2) 


I 


1 


where 


fR  =  f2^,  yk,  zk)  +  (zj  +  z?)  cos («k)  cos  (f^)  yk> 


(6.3) 


-  2(^)  [f(xk,  yk,  zk)  t  Zj  cos(c^.)  cos(^k)]  . 
[sinfP^)  cos(y)  +  cos(/3k)  sinfa^)  sin(y)] 


(6.4) 


and  where; 
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f(Xk»  yk»  zk)  =  ^  sinfc^)  cos(^k)  +  yk  sin(^k)  -  zfc  cos^)  cos(^k).  (6.  5) 


2 

Thus  s  is  a  periodic  function  of  l/r  with  frequencies  fk>  k  =  1, . . . 

-  2 

n,  and  phases  <//k>  k  =  1, . . .  n.  Thus  if  s  is  spectrum  analyzed,  fk  and  i^k  can 
be  determined.  It  can  be  shown  that  by  suitably  varying  y  and  the  para¬ 
meters  ctj^,  ^k  and  f yk>  z.fc)  can  be  obtained.  The  three  parameters  uni¬ 
quely  define  the  kth  specular  multipath  source. 

In  Section  3,  Fig.  19,  we  plotted  a  return  with  n  =  2  as  a  function  of 
l/r.  The  multiple  sinusoids  were  quite  evident.  We  have  plotted  in  Fig.  54 
the  real  part  of  the  Fourier  transform,  in  Fig.  55  the  imaginary  part  and  in 
Fig.  56  the  power  spectrum  of  the  range  corrected  waveform  of  Fig.  19.  The 
three  peaks  in  the  power  spectrum  are  quite  evident.  The  lower  peak  is  due 
to  the  ground  reflection  and  the  upper  peak  is  due  to  the  tilted  plane  reflection. 
The  central  peak  is  a  result  of  beating  between  these  two  multipath  signals. 


6.  3  ANTENNA  DESIGN 

The  effects  of  multipath  strongly  impact  on  the  choice  of  an  antenna  for 
DABS.  The  discussion  of  azimuth  estimation  performance  in  Section  3  demon¬ 
strated  that  the  effect  of  loss  of  SNR  could  be  quite  serious.  This  loss  was  due  to 
the  coherent  interference  between  the  direct  path  and  the  vertical  lobing  path. 

This  effect  can  be  lessened  by  using  an  antenna  with  vertical  aperture  which 
in  effect  will  have  a  cutoff  below  the  horizon.  That  is,  the  signal  coming  in 
at  negative  elevation  angles  will  be  attenuated.  This  results  in  an  effective 
decrease  in  the  reflection  coefficient  and  increase  in  minimum  SNR. 
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REAL  SPECTRUM 


Fig.  54.  Real  par*  ">f  Fourier  transform. 


FREQUENCY 


Fig.  56.  Power  spectrum. 
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Other  effects  can  also  be  studied.  For  example,  the  choice  of  hori¬ 
zontal  aperture  and  the  nature  of  the  sidelobes  is  also  influenced  by  the  type 
of  multipath  present. 

6.  4  SIGNAL  PROCESSING 

The  fade  margins  can  be  determined  by  the  depths  of  the  nulls  that 
occur.  If  these  nulls  are  quite  large  then  more  sophisticated  methods  of  cod¬ 
ing  and  signal  processing  techniques  will  be  necessary. 

6. 5  SITING 

An  important  issue  in  the  determination  of  the  implementation  of  a 
DABS  sensor  is  its  siting.  For  example,  where  should  it  be,  how  high,  and 
what  type  of  tilts  should  be  used.  The  effect  of  specular  multipath  on  these 
issues  is  significant.  Signal  fading  and  coverage  as  well  as  angle  estimation 
performance  are  often  dominated  by  the  multipath  environment.  At  the  pre¬ 
sent,  there  seems  to  be  no  cohesive  measure  of  optimum  siting  performance 
but  possibly,  through  a  detailed  study  of  the  many  multipath  effects,  one  may 
evolve. 


6. 6  CONCLUSIONS 

The  models  that  we  have  proposed  and  d-veloped  can  be  used  to  eval¬ 
uate  and  analyze  the  effects  of  multipath  on  system  performance  and  design. 
Specifically,  as  we  have  observed,  the  effects  of  specular  and  diffuse  multi- 
path  on  azimuth  errors  can  be  evaluated  in  detail.  Furthermore,  link  relia¬ 
bility  can  be  deteriorated  by  specular  fading.  The  effects  of  antenna  gains 
and  cutoffs  in  eliminating  deep  fades  and  serious  azimuth  errors  can  be  ascer¬ 


tained. 


The  localization  algorithm  has  been  found  to  provide  a  vseful  tool 
in  determining  the  source  of  specular  returns.  By  the  nature  of  the  signal 
processing  involved  an  effective  increase  in  aperture  can  be  obtained  thus 
increasing  resolution.  The  use  of  this  technique  on  the  analysis  of  DABSEF 
data  is  presently  under  way. 
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SECTION  SEVEN 


-CONCLUSIONS 


The  ability  to  obtain  azimuth  information  and  communicate  data  over 
an  air-to-ground  surveillance  link  depends  quite  strongly  on  the  different 
types  of  multipath  effects.  To  determine  the  influence  of  multipath,  we  have 
in  this  report,  developed  models  for  the  signals  that  are  received  after  having 
propagated  over  multipath  channels.  The  models  include  the  effects  of  the 
topographical  and  electrical  properties  of  the  reflecting  surfaces  as  well  as 
the  position  of  aircraft  and  sensor. 

As  a  result  of  this  study  we  have  reached  the  following  set  of  conclu¬ 
sions: 

1.  The  most  serious  cause  of  azimuth  angle  of  arrival  estimate 
errors  at  the  ground  terminal  will  be  specular  multipath  sig¬ 
nals  in  the  main  beam.  These  will  tend  to  cause  bias  errors 
in  the  estimate.  The  seriousness  of  the  errors  as  measured 
by  their  magnitude  will  depend  upon  the  angle  made  by  the  ray 
from  the  aircraft  to  the  multipath  reflecting  surface  and  the 
beam  pattern  of  the  receiving  antenna.  If  this  angle  is  small, 
then  the  Fresnel  reflection  coefficient  may  be  near  unity  which 
can  only  be  reduced  by  the  antenna  pattern.  The  magnitude  of 
the  azimuth  estimation  error  also  depends  on  the  total  phase 
difference  between  the  direct  path  and  the  reflected  path.  This 


2. 


3. 


can  be  obtained  from  the  geometry  of  the  aircraft,  sensor  and 
reflecting  surface. 

The  time  delays  for  main  beam  multipath  may  vary  signi¬ 
ficantly,  from  a  few  nanoseconds  to  several  microseconds,  as 
observed  from  DABSEF  data.  Again,  this  depends  upon  geome¬ 
trical  considerations.  The  nature  of  the  surface  will  also  affec* 
the  reflected  signal.  By  using  the  concept  of  Fresnel  zones,  dif¬ 
fraction  effects  can  be  included  in  a  first  order  manner. 

Diffuse  multipath  will  have  a  minimal  effect  upon  ground  ter¬ 
minal  direction  finding  accuracy.  The  analysis  of  Section  4 
has  shown  that  diffuse  multipath  azimuth  errors,  as  measured 
in  standard  deviations,  have  magnitudes  less  than  noise  errors. 
Furthermore,  it  has  been  shown  by  McGarty  [1]  that  diffuse 
multipath  in  certain  cases  contains  information  on  the  azimuth 
of  the  target.  This  is  in  sharp  contrast  to  the  analysis  of 
Kulke,  et  al.  ,  who  employed  results  from  Barton  and  Ward  to 
compute  the  deleterious  effect  of  diffuse  multipath. 

For  most  cases  of  interest,  neither  specular  nor  diffuse  multi- 
path  will  seriously  limit  either  the  direction  finding  or  communi¬ 
cation  capabilities  of  DABS.  There  are  certain  pathological 
situations  concerning  specular  multipath  however,  in  which  the 
effect  will  be  quite  deleterious.  These  have  been  analyzed  in 
detail  by  McAulay. 
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The  second  order  effects,  such  as  shadowing  and  multiple 
oriented  surfaces,  could  possibly  dominate  the  signal  behavior 
at  low  angles.  To  answer  this  question  adequately,  extensive 
experimental  data  are  necessary. 

As  a  result  of  the  preceeding  analysis  it  becomes  evident  that  a  sub¬ 
stantial  experimental  program  should  be  initiated  to  provide  a  data  base 
against  which  the  proposed  models  can  be  compared.  Part  of  this  program 
is  included  in  the  DABSEF  (DABS  Experimental  Facility)  presently  in  opera¬ 
tion  at  Lincoln  Laboratory.  A  more  extensive  program  is  also  envisioned 
which  would  include  a  portable  system. 


APPENDIX  I 


In  the  Appendix,  we  evaluate  in  closed  form  the  bounding  ellipse  dis¬ 
cussed  in  Section  4.  Further,  we  relate  this  to  several  other  physical  para¬ 
meters,  namely,  the  minimum  distance  between  the  source  and  the  receiver 
via  a  multipath  reflection.  The  basic  equality  to  evaluate  the  ellipse,  which 
represents  the  locus  of  all  points  on  the  scattering  plane  whose  total  path 
length  from  receiver  to  transmitter  is  p,  is; 

RY  +  Rz=  p  (LI) 

where  Rj  and  R2  are  given  by  (4.  13)  and  (4.  14)  respectively,  and  p  equals  ct. 
Now  rearrangeing  and  squaring,  we  obtain 

R^  =  p2  -  2p  R2  +  R^  .  (I.  2) 

Substituting  the  values  for  these  quantities  we  obtain, 

2  .  2  2  2  ~  I.  .2  ,  2  ,  2 

x  +  y  +  Zj  =  p  ■  2p  yj{r  -  y)  +x  +  z2 

2,2  -  ,2.2  /T 

+  y  +  x  -  2yr  +  r  +  z  (I.  3) 

Cancelling  and  rearranging,  we  obtain 
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Now  the  term  on  the  right  in  the  above  equality  can  be  simplified  by  writing 


it  as 


f,  2  ,  2  2  ,  2.2  2  2. 

J(p  +  z2  “  1  +  r  *  4<P  "  r  ) 

./  2  2  ,  2  2\  .,  2  2, 

-  4ip  r  +  p  zJ  4(p  -  r  ) 


+  (4p2  r)2  -  2(4p2  r)  (2r(c2  +  z2  -  z2  +  r2)) 


+  4r2{p2  +  z2  -  z2  +  r2)2J  /(4(p2  -  r2)) 


.2,2,  2  2  ,  2. 2  , ,  2  2,2  2  2,  2, 

4p  (p  +  z^  -  Zj  +  r  )  -  16p  r  (p  +  z^  -  z^  +  r  ) 


-  16p2(; 


2  2  2  2  4 

zzp  -z2r  -  r 


)j/(4(P2  - 


,  2 r,  2  2  2  ,  2,2  _ , _  2,  ,  2  2  2  2, 

4p  f(p  +  z2  -  Zj  +  r  )  -  2(2r  )  (p  +  z2  -  zi  +  r  ) 


,  4,  .,  2  4  .,2/  2  2  2  2  4\|/,..  2 

+  4r  ]  -  1 6 p  r  -  I6p  lz2  p  '  z2  r  "  r  ]|/(4(p  - 


.  2(2  2  .  2  22  .,2  2,  2  2. 

4p  p  -  r  +  z2  -  Zj  -  1 6 p  z2(p  -  r  ) 

- - - 

4(p*  -  r  ) 


(1.8) 


which  is  the  desired  result. 
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Now  we  can  ask  the  question,  where  is  the  point  at  which  the  mini¬ 
mum  time  ray  touches  the  surface?  That  is  for  a  given  p,  the  time  from 
source  to  receiver  is  (R^  -  R^/c.  We  first  want  to  find  the  value  of  x  which 
minimizes  this  expression.  By  standard  differentiation  techniques  we  obtain 


x  . 
min 


r  z, 


z1+  z2 


(1.9) 


Note  that  if  z,  =  z.,,  then  x  .  =  d/2.  Now  the  time  it  takes  to  cover  this  mini 

i  z  min 

mum  distance  is  t^,  where 


(I.  10) 


Compare  this  to  the  time  between  receiver  and  transmitter,  t^, 

t  *“  ' 

r  c 


(I.  11) 


Now  at  t  =  tg,  p  equals  c  tQ. 

Then  substitution  into  vQ  yields  rQ  equal  to  zero  and  xQ  equal  to  xmlQ. 
Thus  as  t  increases  beyond  tQ,  the  ellipse  forms  and  grows. 
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antenna  area 

area  of  first  Fresnel  zone 
area  of  Lth  reflector 
irradiated  area 

complex  envelope  of  time  invariant  portion  of  diffuse  field 

complex  envelope  of  diffuse  field 

velocity  of  propagation 
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divergence  factor 

expectation 
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current  density 
monopulse  slope 
plane  wave  vector 

components  of  propagation  vector 

channel  spread  function!  azimuth  only 
channel  spread  function 
covariance  of  diffuse  signal 
array  delay  vector 
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number  of  array  elements 
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range 
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scattering  surface 
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interference -to -signal  ratio 
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channel  attenuation  coefficient 
standard  deviation  of  surface  height 

diffuse  signal  energy- 

power  in  diffuse  field 
s catte ring  cross-se ction 
conductivity  of  ith  scattering  plane 
correlation  length  of  surface 

noise  covariance 

signal  energy 

scattering  cross  sections 
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variance  of  azimuth  estimate 


2 

a 

H- 


T 


4>i 

<PF 


%i 


variance  of  p 

multipath-interference  statistic 
elevation  (incident  or  ith  signal) 
elevation  of  scattered  signal 
phase  due  to  Fresnel  reflection 

phase  difference 

multipath  phase 
phase  due  to  range 

carrier  frequency- 
array  ambiguity  function 
transformation  matrix 


144 


REFERENCES 


Ament,  W.S. 

[1]  "Toward  a  Theory  of  Reflection  by  a  Rough  Surface,  "  Proc.  IRE, 

4R  (1953)  pp.  142-146. 

[2]  "Forward  and  Back-Scattering  From  Certain  Rough  Surfaces,  " 

KE,  AP-4,  (1965)  pp.  369-373. 

Barrick,  D.  E. 

[1]  "Rough  Surface  Scattering  Based  on  the  Specular  Point  Theory,  " 
IEEE,  AP-16,  No.  4,  (July  1968)  pp.  449-454. 

[2]  "Unacceptable  Weight  Correlation  Coefficients  and  the  Quasi- 
Specular  Component  in  Rough  Surface  Scattering,  "  Radio  Science 
Vol.  5,  No.  4,  (April  1970)  pp.  647-654. 

[3]  "Theory  of  HF  and  UHF  Propagation  Across  the  Rough  Sea,  " 

Radio  Science,  Vol.  6,  No.  5,  (May  1971). 

Barrick  D.  E.  ,  Peake,  W.  H. 

[1]  "A  Review  of  Scattering  From  Surfaces  with  Different  Roughness 
Scales,  "  Radio  Science,  Vol.  3,  No.  8,  (August  1968)  pp.  865-868. 

Barton,  D.  K,  ,  Ward,  H.  R. 

fl  1  Handbook  of  Radar  Measurement,  (Prentice -Hall,  Englewood  Cliffs, 

to;;  'TT6'9y; - 

Bass,  F.  G.  ,  et  al. 

[1]  "Very  High  Frequency  Radio  Wave  Scattering  by  a  Distributed  Sea 
Surface  Pt.  I,  "IEEE,  Vol.  AP-16,  No.  5,  (Sept.  1968)  pp.  554- 
559. 

[2]  "Very  High  Frequency  Radio  Wave  Scattering  by  a  Distributed  Sea 
Surface  Pt.  II,  "  IEEE,  Vol.  AP-16,  No.  5  pp.  560-588. 

Beard,  C,  I. 

[1]  "Coherent  and  Incoherent  Scattering  of  Microwaves  from  the  Ocean 
IEEE,  Vox.  AP-9,  No.  5,  (Sept.  1961)  pp.  470-483. 

[2]  "Behavior  of  Non-Rayleigh  Statistics  of  Microwave  Forward  Scatter 
from  a  Random  Water  Surface,  11  IEEE,  Vol.  AP-15,  No.  5 

(Sept.  1967)  pp.  649-657. 


145 


Beard,  C.I.  ,  etal. 

[1]  "Phenomenological  Vector  Model  of  Microwave  Reflection  from  the 
Ocean,  "IRE,  Vol.  AP-4,  (1956)  pp.  162-167. 

Beckmann,  P. 

[1]  "Scattering  by  Composite  Rough  Surfaces,  "  Proc.  IEEE, 

(Aug.  1965)  pp.  1012-1015. 

[2]  "Scattering  of  Light  by  Rough  Surfaces,  "  Progress  in  Optics, 

IV,  (E.  Wolf,  Ed.,  J.  Wiley,  N.  Y.  ,  1967)  pp.  55-70. 

[3]  "Shadowing  of  Random  Rough  Surfaces,  "  IEEE,  AP-14,  (1965) 
pp.  384-388. 

Beckmann,  P.  ,  Spizzichino,  A. 

[1  ]  The  Scattering  of  Electromagnetic  Waves  from  Rough  Surfaces, 
(Permagon,  N.  Y.  ,  19^3^ 

Brockelmann,  R.  A.  ,  Hagfors,  J. 

[1]  "Note  on  the  Effect  of  Shadowing  on  the  Backscattering  of  Waves 
from  a  Random  F.ough  Surface,  "  IEEE,  Vol.  AP-14,  No.  5, 

(Sept.  1966)  pp.  621-629. 

Clarke,  R.  H.  ,  Hendry,  G.  O. 

[1]  "Prediction  and  Measurement  of  the  Coherent  and  Incoherent 
Power  Reflected  from  a  Rough  Surface,  "  IEEE,  Vol.  AP-13, 

No.  3,  (1964)  pp.  353-363. 

DeLorenzo,  J.  D.  ,  Cassedy,  E.S. 

[1]  "A  Study  of  the  Mechanism  of  Sea  Surface  Scattering,  "  IEEE, 

Vol.  AP-14,  No.  15,  (Sept.  1966)  pp.  611-620. 

Eckart,  C. 

[1  ]  "The  Scattering  of  Sound  from  the  Sea  Surface,  "  J.  A.  S.  A.  , 

Vol.  25,  No.  3  (May  1953)  pp.  566-570. 

Fung,  A.  K.  ,  Moore,  R.  K. 

[1]  "The  Correlation  Function  in  Kirchoff  's  Method  of  Solution  of 

Scattering  of  Waves  from  Statistically  Rough  Surfaces,  "  J.  Geo. 
Res.,  Vol.  71,  No.  12,  (June  1966)  pp.  29  3  9  -  29  43. 

Furutsu,  Kc 

[1]  "A  Statistical  Theory  of  Range  Diffraction,  "  Radio  Sci.  ,  Vol.  1, 

No.  1,  (Jan.  1966)  pp.  79-98. 


146 


Hagfors,  T. 

[1]  "Backscattering  from  an  Undulating  Surface  with  Applications  to 
Radar  Return  from  the  Moon,  "  J.  Geo.  Res.,  Vol.  69,  No.  18, 
(Sept.  1964)  pp.  3779-3784. 

Hayre,  H.  S. ,  Moore,  R.  K. 

[1]  "Theoretical  Scattering  Coefficient  for  Near  Vertical  Incidence 
from  Contour  Maps,  "  Jour.  Res.  Nat.  Bur.  Sts.  ,  Series  D, 

Vol.  65D,  No.  5,  (Sept. -Oct.  1961)  pp.  427-432. 

Hoffman,  W.C. 

[1]  "Scattering  of  Electromagnetic  Waves  from  a  Random  Surface,  " 
Quant.  Appl.  Math.,  J_3,  (1955)  pp.  291-304. 

Jordan,  E.  C. ,  Balmain,  K.  G. 

[1]  Electromagnetic  Waves  and  Radiating  Systems,  (Prentice  Hall, 
ITnglewood  Cliffs"!  TO!  7  T«568 j: - - - 

Kerr,  U.  E. 

[1]  Propagation  of  Short  Radio  Waves,  (Dover,  N.  Y.  ,  1965). 

Kodis,  R.  D. 

[1]  "A  Note  on  the  Theory  of  Scattering  from  an  Irregular  Surface,  " 
IEEE,  Vol.  AP-14,  No.  1,  (Jan.  1966)  pp.  77-82. 

Lonquet-Higgins,  M.  S. 

[1]  "Reflection  and  Refraction  at  a  Random  Moving  Surface,  PL  I,  " 
J.O.S.A.  ,  Vol.  50.,  No.  9,  (Sept.  I960)  pp.  838-844. 

[2]  "Reflection  and  Refraction  at  a  Random  Moving  Surface,  Pt.  II,  " 
J.O.S.A.,  Vol.  5,  No.  9,  (Sept.  I960)  pp.  845-850. 

Maltz,  F. H.  ,  Meecham,  W.C. 

[1]  "The  Effect  of  Curvature  Variations  on  the  Scattering  from  Rough 
Surfaces,  "J,  Stat.  Phy. ,  Vol.  3,  No.  4,  (1971)  pp.  385-394. 

McAulay,  R.  J. 

[1]  "The  Effects  of  Interference  on  Monopulse  Performance"  Lincoln 
Labortory,  M.  I.  T. ,  Technical  Note,  TN-1973-30,  (1  Aug.  1973) 
ESD-TA-73-176. 


McAulay,  R.  J.  ,  McGarty,  T.  P, 

[1]  "An  Optimum  Interference  Detector  for  DABS  Monopulse  Data 

Editing,  "  Lincoln  Laboratory  M.  I.  T.  Technical  Note,  TN-1973-48 
(26  Sept.  1973)  ESD-TR-73-253. 


McGarty,  T.  P. 

[1]  "Azimuth-Elevation  Estimation  Performance  of  A  Spatially 

Dispersive  Channel,  "  IEEE  Vol.  AES-10,  No.  1  (January  1974) 
pp.  58-69. 

Mitzner,  K.  M 

[1]  "Change  in  Polarization  on  Reflection  from  a  Tilted  Plane,  "  Rad. 

Sci.  Vol.  1,  No.  1,  (Jan.  1966)  pp.  27-29. 

Rice,  S.  O. 

[1]  "Reflection  of  Electromagnetic  Waves  from  Slightly  Rough  Surfaces,  " 
Comm.  Pure  &  Appl.  Math,  4,  (1951)  pp.  351-378. 

Silver,  S. 

[1]  Microwave  Antenna  Theory  and  Design,  (Dover,  N.  Y.  ,  1965). 

Smith,  B.  G. 

[1]  "Geometrical  Shadowing  of  a  Random  Rough  Surface,  "  IEEE, 

Vol.  AP-15,  No.  3,  (Sept.  1967)  pp.  668-671. 


Spingler,  G.  F. 

[1]  "Analysis  of  Side-Lobe  Suppression  Capabilities  at  Joint-Use  Radar 
Beacon  Sites,"  FAA  Report  No.  RD-67-59. 

Stratton,  J.  A. 

[1]  Electromagnetic  Theory,  (McGraw  Hill,  N.  Y. ,  1941). 

Twersky,  V. 

[1]  "Reflection  Coefficients  for  Certain  Rough  Surfaces,  "  J.  App.  Phys.  , 

24,  (1953)  pp.  659-660. 

[2]  "Certain  Transmission  and  Reflection  Theorems,  "  J.  App.  Phys. , 
Vol.  25,  No.  7,  (July  1954)  pp.  859-862. 

[3]  "On  the  Scattering  and  Reflection  of  Sound  by  Rough  Surfaces,  " 
J.A.S.A.  Vol.  29,  No.  2,  (Feb.  1957)  pp.  209-225. 

[4]  "On  Scattering  and  Reflection  of  Electromagnetic  Waves  from  Rough 
Surfaces,  "  IRE,  AP-5,  (1957)  pp.  81-90. 

Van  Trees,  H.  L. 

f  1 7  Detection  Estimation  and  Modulation  Theory,  Pt.  m,  (J.  Wiley,  N.Y. 

T9TTTT - 

Venetsanopoulos,  A.  N. ,  Tuteur,  F.  B. 

[1]  "Stochastic -Filter  Modeling  for  the  Sea-Surface  Scattering  Channel,  " 
J.A.S.A. ,  49,  (1971)  pp.  1100-11007. 


V 


Wagner,  R.  J. 

[1]  "Shadowing  of  Randomly  Rough  Surfaces,  "  J.A.S.A.,  Vol.  41, 
No.  1,  (1967)  pp.  138-147. 

Yaglom,  A.M. 

[1]  "Second-order  Homogneous  Random  Fields,  "  Proc.  4th  Berkeley 
Symp.  Math.  Stat.  and  Prob.  2,  pp.  593-620. 


149 


*  U.S.  GOVERNMENT  PfUNTONG  OFFICE:  J9T4  1IM-1M  SO 


